Analysis of phonon-induced spin relaxation processes in silicon 
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We study all of the leading-order contributions to spin relaxation of conduction electrons in silicon 
due to the electron-phonon interaction. Using group theory, k -p perturbation method and rigid- ion 
model, we derive an extensive set of matrix element expressions for all of the important spin-flip 
transitions in the multi- valley conduction band. The scattering angle has an explicit dependence on 
the electron wavevectors, phonon polarization, valley position and spin orientation of the electron. 
Comparison of the derived analytical expressions with results of empirical pseudopotential and 
adiabatic band charge models shows excellent agreement. 
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I. INTRODUCTION 

Silicon is an ideal material choice for semiconductor 
spin-based devices^ - — It has a relatively weak spin-orbit 
coupling which leads to a negligible probability of flipping 
the electron's spin during a scattering event. Further- 
more, the Dyakonov-Perel spin relaxation mechanism is 
absent in bulk silicon due to its inversion symmetry (i.e., 
no intrinsic magnetic field around which the electron spin 
precesses) A Finally, the zero nuclear spin of its naturally 
abundant isotope suppresses spin relaxation by hyperfine 
interactions.— ~— These characteristics have motivated a 
wide interest in silicon spintronics£~— 

In cases where spin dephasing by precession is weak, 
the spin relaxation rate is set by spin-flip scattering. 
For conduction electrons in crystals, this scattering is 
described by a matrix element of the form (k.2 > JJ-n 
\H S { |ki, fl n ) where the states are identified with wavevec- 
tor k due to the translational symmetry For a given 
spin- flip mechanism (H s {), this matrix element depends 
on the initial and final state wavevectors ki and k2, as 
well as on the spin orientation n. In silicon, conduc- 
tion electrons reside in six valleys near the edges of the 
Brillouin zone as shown in Fig. [IJa). The valley cen- 
ters are positioned on the A-axis connecting the T and 
X points. In multivallcy semiconductors, the electron 
remains in the same valley after an intravalley scatter- 
ing and it switches valleys after an intervalley scattering. 
Figure [TJ(a) shows examples of intervalley scattering in 
Si, where g and / processes refer to electron scattering 
between opposite valleys and between valleys that re- 
side on perpendicular crystal axes, respectively. When 
dealing with spin-flip scattering, the spin orientation (n) 
further breaks the scattering symmetry. The resulting 
anisotropy in spin relaxation depends on the projections 
of n on crystallographic axes [Fig. QJb)]. 

Spin-flip mechanisms can be classified into Yafet and 
Elliott processes. The former involves spin-dependent 
interaction whereas the states are viewed as pure spin 
states.— Examples of spin-dependent interactions include 
the spin-orbit coupling of the host crystal via electron- 
phonon interaction (H s f ~ 5R ■ Vy so ), spin-orbit cou- 




FIG. 1: (Color online) (a) Valley positions and high symmetry 
points in the Brillouin zone. Valley centers are about 0.15 x 
2n/a from the X points where a is the lattice constant. Also 
marked are representative intervalley transitions by g and / 
processes, (b) The spin orientation (n) with respect to the 
crystallographic axes. 



pling of defects via electron-impurity scattering (H s { ~ 
VVi m p x p • a), and electron- nuclear hyperfine interac- 
tion (H s { ~ As • I) . The Elliott processes are governed by 
spin mixing in the electron states due to the crystal spin- 
orbit coupling, whereas the interaction is spin indepen- 
dent (e.g., SR ■ W), 1 - In this paper the focus is on spin 
relaxation due to the electron-phonon interaction where 
the interplay between Elliott and Yafet processes plays a 
key role in setting the intrinsic spin lifetime. In Sec.HIlwe 
provide a general overview of the spin relaxation mech- 
anisms in n-type silicon and we show that the electron- 
phonon interaction dominates the spin relaxation of con- 
duction electrons over a wide range of temperature and 
doping conditions. 

The spin relaxation due to electron-phonon interac- 
tion in silicon was studied using different theoretical 
approaches. The seminal works of Ovcrhauser, Elliott 
and Yafet enabled an approximate quantitative connec- 
tion between the spin and momentum relaxation times 
via the shift of the g-i&ctor Yafet also derived a 
general form of the intravalley scattering by acoustic 
phonon modest Intervalley scattering was included in 
three recent works. 20-22 Numerical analysis of the mca- 
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sured spin relaxation times above 50 K in intrinsic silicon 
was carried out by Cheng, Wu and Fabian^ The anal- 
ysis included calculation of electron energies and states 
using an empirical pseudopotential model (EPM)^ of 
phonon dispersion and polarization vectors using an adi- 
abatic bond charge model (ABCM),— and of a rigid- 
ion model approximation to describe the electron-phonon 
interaction^ An analytical approach was then devel- 
oped by Li and Dcry using a compact spin-dependent 
k-p Hamiltonian model around the zone edge X point^ 
This model was then used to derive the dominant spin- 
flip matrix elements where the spin-orbit coupling signa- 
ture appeared only in the expansion of electronic states. 
The electron-phonon interaction, on the other hand, was 
mimicked by a phenomenological connection with defor- 
mation potential and scattering parameters. In another 
recent theory, Tang, Collins and Flatte have used a tight- 
binding model to calculate the spin relaxation in strained 
silicon and germanium^ 

In this paper we present a comprehensive theory of 
phonon-induced spin-flip mechanisms of conduction elec- 
trons in bulk silicon. Matrix elements are derived us- 
ing two different approaches based on whether they are 
wavevector dependent or not. TableUpresents the power- 
law dependence for spin-flip and spin-conserving matrix 
elements as well as the theoretical approaches to de- 
rive their forms. When transitions between valley-center 
states do not vanish [i.e., M(0,si ; 0,s 2 ) + 1 0], group the- 
ory can be used alone to derive explicit forms of their 
matrix elements including the spin orientation depen- 
dence. As shown in the right column of the table, spin- 
conserving scattering is wavevector independent other 
than with long-wavelength acoustic phonons. On the 
other hand, spin-flip scattering is wavevector indepen- 
dent only for the /-process whereas other processes have 
higher power-law dependence than in spin-conserving 
scattering. Accordingly, the /-process is dominant in 
spin relaxation whereas its weight is comparable to other 
processes in momentum relaxation. In the next three 
paragraphs we discuss key aspects of the present work in- 
cluding explanation of the approach to derive wavevector- 
dependent matrix elements. 

For every type of electron-phonon scattering, we ex- 
press the leading-order spin-flip matrix elements as func- 
tions of the electron wavevector (k), phonon polariza- 
tion (£), valley position of the electron and its spin ori- 
entation (n). For the important /-process we derive a 
complete set of selection rules by rendering double group 
representation matrices in conjunction with time reversal 
symmetry. This approach is more informative than the 
common practice of using the character table since it dis- 
tinguishes between momentum and spin scattering pro- 
cesses, it unambiguously identifies parameters for both 
processes, and it is generalized for any spin orientation. 
The latter is important in analyzing experimental mea- 
surements where the orientation of injected spins is dic- 
tated by the shape and magnctocrystallinc anisotropy of 
ferromagnetic contacts or by the propagation and hclicity 



TABLE I: Power-law dependence of leading order matrix ele- 
ments for all types of electron-phonon scattering. ki,2 are the 
initial and final electron wavevectors and they are measured 
from their respective valley centers, q and K denote their 
difference and average, respectively. Intravalley scattering is 
divided to interaction with long-wavelength optical (OP) and 
acoustic (AC) phonons. Also mentioned are theoretical ap- 
proaches to derive explicit forms of the matrix elements (be- 
yond the power-law dependence). See text for further details. 





M(ki,s; k 2 ,-s) 


M(ki,s; k 2 ,s) 


/-process 


k- independent, 
double group. 


k-independent, 
single group. 


(/-process 


linear in K, 
'k-p'+single group. 


k-independent, 
single group. 


intravalley 


AC: quadratic in q, 
OP: linear in q, 
'k-p'+single group. 


AC: linear in q, 
OP: k-independent, 
'k-p'+single group. 



of a circularly-polarized light beam. 

In the analysis of intravalley and g-process spin flips 
we employ a different approach. Derivation of selection 
rules in these cases cannot be carried out solely by group 
theory because the spin-flip matrix elements are wavevec- 
tor dependent. Instead, we utilize a combined approach 
that involves spin-dependent k ■ p, rigid-ion and group 
theories. Ad hoc selection rules are derived for electron- 
phonon interaction between basis functions that appear 
in the k-p expansion of states around the valley center. 
The selection rules are derived from single group theory 
whenever multiple bands have to be included in the k • p 
expansion. Double group bares no advantage in this case. 
This theoretical procedure allows us to include the crys- 
tal symmetry not only in the electronic states but also in 
the electron-phonon interactions. One outcome is that 
new scattering-angle symmetries are revealed. We derive 
appealing matrix element forms by employing an elas- 
tic continuum approximation for diamond crystal struc- 
tures. A second outcome is that the out-of-phase motion 
of atoms in the primitive cell is quantified and shown to 
play a role in all types of intravalley spin-flip processes 
including scattering with acoustic phonons. 

In this paper we thoroughly study the proximity effect 
between the lowest pair of conduction bands in silicon. 
In his seminal work, Yafet showed that the intravalley 
spin-flip matrix element due to scattering with acous- 
tic phonon modes is proportional to the general form 
C A so q / Eg where C is a deformation potential constant, 
Eg is the energy gap from bands where the spin-orbit cou- 
pling is considered (upper valence bands), and A so is the 
strength of the coupling^ The quadratic dependence on 
the phonon wavevector (q 2 ) is a consequence of time re- 
versal and space inversion symmetries. This quadratic 
dependence allowed Yafet to predict the T -5 / 2 temper- 
ature dependence of the spin lifetime due to intravalley 
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electron scattering with acoustic phonon modes. Yafct's 
general matrix element, however, misses two important 
features. First, the scalar form of this matrix element 
does not provide details of the dominant phonon modes 
and cannot capture the dependence on directions of q 
and n. Second, the interband deformation potential be- 
tween the lowest pair of conduction bands plays a crucial 
role in intravalley spin relaxation of silicon. As a result, 
Ago/Eg should be replaced by A so /(E g Ac) where Ac 
denotes the band gap between the two conduction bands 
at the valley center^ This interband coupling is a unique 
feature of the conduction band in silicon. It results from 
the proximity of the valley center to the two-band de- 
generacy at the X point. This coupling also explains the 
discrepancy between theoretical predictions of the spin 
lifetimes in strained siliconi 15 ' 22 This discrepancy can be 
resolved by inclusion of d-orbitals in the tight-binding 
model to correctly capture the electronic behavior near 
the X points 



Paper organization and outline of central results 

This paper is organized as follows. Section [TTI surveys 
various spin relaxation processes in n-type silicon. Sec- 
tion |TlT] provides a theoretical framework for the electron- 
phonon interaction. The dominant /-process is stud- 
ied in Sec. IIV1 Using four scattering constants, Eq. (fl7|) 
lists explicit forms of the /-process matrix elements as 
a function of the phonon symmetry and of the spin ori- 
entation. In section [V] we present detailed derivations 
and extensive results of intravalley and <?-process spin 
flips. This section begins with an introduction of the 
special features and tabulation of representative results 
(Table |V| . Selection rules and Hamiltonian forms are in- 
troduced in section IV Al as well as accurate state vec- 
tors up to quadratic-in-k terms. Sections IV BIIVDl in- 
clude derivations of the core spin-flip matrix element ex- 
pressions with spin orientation dependence [Eqs. (|57[1 - 
(|59| with their deformation potential constants defined 
in Eqs. (|ig]). (gSJ and ([55])]. Some insights about the 
relation between analytical derivation and EPM are dis- 
cussed in Scc lVIl Spin lifetimes from integrating analyt- 
ical matrix elements [Eqs. (|oTj) - (pKt]) ] as well as different 
levels of numerical calculations are presented and com- 
pared in Section IVlII Section lVIIII includes a summary of 
results and an outlook on future work. Appendices lAlFl 
include technical details from different derivation phases. 



II. SPIN RELAXATION PROCESSES 

Spin lifetime of electrons in n-type bulk silicon varies 
by more than 12 orders of magnitude with changing tem- 
perature and doping conditions.— ~— We summarize the 
important mechanisms using the map in Fig. [2] 

In the low-temperature and non-degenerate regime (re- 
gion 1), electrons freeze on isolated donor sites. The 
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FIG. 2: (Color online) Diagram of dominant spin relaxation 
mechanisms in n-type silicon as a function of temperature and 
donor concentration. It is applicable at or near equilibrium 
conditions. In region 1 electrons are localized on isolated im- 
purity sites. In regions 4(3) they populate the conduction 
(impurity) band. Region 2 is a precursor of the impurity 
band (donor clusters). Region 5 includes more than a sin- 
gle phase. Transitions between regions are generally gradual 
and colors as well as dashed border lines are added for ease 
of illustration. The findings in this paper are relevant for 
non-degenerate conditions (<10 17 cm -3 ) in region 4. 



ground states is split into 6 states (without spin) due to 
valley-orbit coupling.— For typical shallow donors such 
as phosphorus the singlet (nondegenerate) state is lo- 
cated ~45 meV below the conduction band edge, while 
the doublet (two-fold degenerate) and triplet (three-fold 
degenerate) states are only slightly split and located 
~35 meV below the conduction band edge. The ex- 
tremely long spin lifetime of localized electrons in region 1 
is governed by electron-phonon Raman processes and by 
hyperfine interactions with the non-zero nuclear spin of 
the impurity or Si 29 isotopes ! 33 ' 50 " — As an example, we 
mention the Orbach process which was shown by Castner 
to dominate between 5 K and 20 K.— Spin flips at the 
low-energy singlet state are caused by phonon-induccd 
virtual transitions to intermediate triplet states at which 
the spin-orbit coupling admixes spin-up and spin-down 
components. This process requires an initial absorption 
of a ^10 meV phonon to mediate a singlet-to-triplet tran- 
sition followed by phonon emission to transfer the elec- 
tron back to the singlet state but with opposite spin. 

The spin relaxation in region 3 of Fig. [5] is due to the 
formation of an impurity band. At these low temper- 
atures and intermediate donor concentrations, the im- 
purity band is populated and separated from the nearly 
empty conduction band. Compared with region 1, the 
spin relaxation is enhanced by orders of magnitude due 
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to the overlap of wave functions in different impurity 
centers^ On the insulator side (region 2), the spin re- 
laxation is governed by exchange interactions between 
localized electrons and by rich hyperfine configurations 
of donor clusters! 29 ' 35 ' 54 In barely metallic samples (re- 
gion 3) the spin admixture is increased by broadening of 
the triplet and singlet bands i 47 ' 55 Spin relaxation is then 
governed by exchange and motional narrowing when elec- 
trons travel across the (random) potential] 39 ' 43 ~ 47 ' 55 

Spin relaxation in region 4 of Fig. [2] is of conduc- 
tion electrons. Transition into region 4 is made cither 
by increasing the temperature (electrons are thermally 
excited to the conduction band), or when the impurity 
band is merged into the conduction band (heavily de- 
generate doping). Region 4 is pertinent to spintronic 
devices in which electrons are swept away from the mag- 
netic junction . 48 ' 5 6 ~ — Then, the relaxation is governed by 
impurity scattering (degenerate doping) or by electron- 
phonon interactions (non-degenerate doping). In lateral 
spintronic devices, on the other hand, electrons are kept 
at the vicinity of the contact ^r— and the transport is af- 
fected by the geometry and properties of the contact i 7 - 3 - - — 
Spin lifetimes that exceed the 1 ns scale were recently re- 
ported in lateral devices that incorporate a highly degen- 
erate n-type silicon channel*^— These reports suggest 
that impurity scattering in the bulk silicon channel plays 
an important role in heavily degenerate lateral devices 
that rely on spin accumulation. 

Electron-phonon interaction plays a key role in the 
non-degenerate doping regime of region 4. In nearly in- 
trinsic conditions, it dominates the spin relaxation al- 
ready above ~30 while at ^10 17 cm~ 3 it dominates 
above ~100 K. The increase in temperature is brought by 
the extended range at which freeze-out conditions persist, 
by exchange between free and localized electrons, and 
by scattering from ionized impurities.— For transport at 
large applied electrical fields (where Fig. [2] is no longer 
valid) , the electron-phonon scattering dominates the spin 
relaxation of conduction electrons at all lattice tempera- 
tures since the effective temperature of drifting electrons 
is significantly higher than the lattice temperature! 87 ' 88 

Electron-impurity scattering dominates the spin relax- 
ation of conduction electrons in highly degenerate dop- 
ing concentrations (top part of region 4). An impor- 
tant aspect is that at intermediate and degenerate con- 
centrations the impurity type influences the spin life- 
time (irrespective of the temperature) i 42 ' 43 ' 55 The origi- 
nal Elliott- Yafct theory, on the other hand, suggests that 
the spin-orbit coupling of the host crystal (rather than 
of the impurity) determines the spin relaxation of con- 
duction electrons! 17 ' 18 Finally, in region 5 several pro- 
cesses coexist and provide measurable spin relaxation. 
At donor concentrations of I0 15 - 10 17 cm~ 3 , Lepine 
showed that the spin relaxation is governed by exchange 
between free and localized electrons, by the larger spin- 
orbit coupling of triplet states, and by the modulation of 
the hyperfine field (when electrons make thermal transi- 
tions between localized states)^ At lower donor concen- 



trations, the role of electron-phonon scattering increases 
since the relative fraction of electrons in the conduction 
band increases ! 38 ' 89 



III. ELECTRON-PHONON INTERACTION 

Within the harmonic approximation framework 
the electron-phonon interaction is generally described 

H p = ]T5R JQ -VR 3Q V at (r-R JQ ), (1) 

where j sums over N primitive cells and a — {A, B} 
labels the two atoms of a primitive cell with the origin 
chosen at their mid-point position, r = ta = —tb = 
— (a/8)(l, 1, 1). An atom position is then denoted by 
Rj a = Rj + T a , and its potential by 

V at (r) = (r)Z + [W a t (r) x p] • a . (2) 

1 is the 2x2 identity matrix and rr is the vector of Pauli 
matrices. The first (second) term is the bare potential 
(spin-orbit coupling), and it corresponds to the Elliott 
(Yafet) part of T-L ep - The atom vibration in Eq. ([T]) is 
expressed as 

6R ja = v / WRqW[aq^(q)e !qE] ° + h-c] , (3) 
q 

where p is the material density and w(q) is phonon fre- 
quency. a q is the annihilation operator, £ is the nor- 
malized phonon polarization vector,— and 'h.c' stands 
for Hermitian conjugate. The square-root prcfactor is 
written in accordance with ^ Q |£a| 2 = 2 (the number 
of atoms in the primitive cell). Using crystal momen- 
tum conservation together with identities of phonon cre- 
ation, (ji2\a '|ni) = \Jn\ + 15 n2 ,n 1 +i, and annihilation, 
(ri2|a|rii) = ^/ni# n2)ni _i, we write the transition ampli- 
tude of an electron from wavevector ki to k2 by, 

(k 2 ,s 2 ;n(q) ± l|"H ep |ki, si; n(q)) = (4) 

~\l M^V n(q) + \ ± I X (^I^pI^Sx) , 
where 

#cp = & (q) eiqRjQ ■ VrVat ( r - R ^), (5) 

j,as 

q = k 2 — ki, n(q) is phonon occupation (Bose-Einstein 
distribution at thermal equilibrium), and s = {f|~,-IJ-} de- 
notes the spin-up and spin-down states. Due to time 
reversal and space inversion symmetries each band at 
wavevector k is spin degenerate and we can define its 
states with respect to the spin orientation n, such that 

(k, n, f \a ■ n| k, fi, f) = -(k, /i, 4 \tr ■ n| k, fj,,if) > , 
(k,/i,ft| CT -n|k )M ,^)=0, (6) 
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TABLE II: In-phase and out-of-phase phonon polarization 
vectors in the long- wavelength limit (q <C 2ir/a). We have 
used the elastic continuum approximation for diamond crys- 
tal structures. A = {TAi,TA 2 ,LA} and Y xyz w a/16. Two 
TA or TO polarizations can be linearly combined into any 
other orthonormal ones. See text for further explanation. 



£tA! > Srch (<l) (Qv,-<l*,0)/y/qZ + q$ 

£ta 2 > 5to 2 (<0 (qxqz,q y qz, -ql - q 2 y )l(\Jqi + <?g|q|) 

^LA.^Lo(q) q/kl 



where /i is the band index. Mixed by spin-orbit interac- 
tion, the defined spin-up and spin-down states (ffj-IJ-) are 
not pure spin states (t,l)- 

Based on the symmetry of electron and phonon states 
we will derive selection rules for M (ki, si; k 2 , s 2 ) = 
(k2,S2|/f e p|ki,si). It is often convenient to work with 
interaction that has a fixed parity under space inversion 
operation of the crystal. We convert Eq. (JT|) into in-phase 
and out-of-phase parts, 

H cp = -^<5R± -V r V± (r - Rj), 

j,± 

SRf = \ (STL, A ± SR jB ) , 
V±(r) = V Bt (r-T. A )±V Bt (r-TB). (7) 

The (anti)symmetrized cell potential V+(_) is associ- 
ated with the in(out-of)-phase motion. Using the crys- 
tal translational symmetry^ £\ e ,<lRj VVi(r-Rj) cou- 
pled between states at ki and k 2 is further reduced to 
A r W±(r)5 q+ k 1 -k 2 .g- With Eq. (0, the matrix element 
becomes 

M(ki,si;k 2 ,s 2 ) = (k 2 , s 2 |-ff ep |ki, Si) 

= iV^^ ± (q)-(k 2)S2 |V r V ± (r)|k 1 ,s 1 ), (8) 

± 

where in-phase and out-of-phase polarization vectors 

€ ± (q) = [eA(q)e*™ ± | B (q)e tqTB ]/2 , (9) 

satisfy |£ + | 2 + |£~| 2 = 1. In what follows we abbreviate 
the notation and absorb the scaling factor into the matrix 
element definition (i.e., N(...) —> (...)). 

In the analysis of intravalley scattering, explicit forms 
of polarization vectors are utilized. To gain further in- 
sight of their symmetries we invoke the elastic continuum 
approximation for diamond crystal structures (incorpo- 
rating internal displacement) Table [TT1 lists forms of 

(q) m t ne long-wavelength limit where A is the phonon 
mode. The in-phase polarization vectors of acoustic 
modes (first three rows) are of the order of unity while 
their out-of-phase vectors (fourth row) are of the order 



of qa/2n <C 1. They flip roles for optical modes (i.e. 
order of unity for out-of-phase and the order of qa/2ir 
for in-phase vectors). Components of out-of-phase vec- 
tors relate to the in-phase components and embody the 
structure of the two-atom primitive cell, 

C(q) = \itM)(^ TA -^' TB ) + *E Fmneqm^M 

Here, A denotes any of the acoustic modes and I, m, n 
are cyclic permutation of the coordinates. In writing 
the expression in the fourth row of the table we have 
used the following approximations. The exponential 
term, brought by the macroscopic strain, is replaced by 
iq • r£ + (q) due to the long- wavelength nature (q — > 0). 
For the other term, brought by internal displacement, 
the only non- vanishing components of the third-rank ten- 
sor Tijk are T xyz = T yzx = T zxy in a diamond-crystal 
structure.*^ Their value is 0.5a/8 . 24 i 97 

It is emphasized that in this work we use the rigid- 
ion model to derive general results of intravalley spin- 
flip matrix elements [Eqs. (|57|) - ([55|) ]. These results do 
not depend on the approximate forms of £ ± in Table [TT1 
However, since these approximations are fairly accurate 
and become exact along high symmetry directions, we 
can make use of them to derive appealing spin-flip matrix 
elements of long- wavelength phonon modes (Table [Vj) . 



IV. INTERVALLEY /-PROCESS SPIN FLIPS 

At elevated temperatures, spin relaxation in un- 
strained silicon is largely dominated by the intcrvalley 
/-process! 15 ' 20 i 22 ' 87 Therefore, it is crucial to have a com- 
plete set of selection rules for /-process spin scatter- 
ing. Selection rules of the intcrvalley momentum relax- 
ation (neglecting spin-orbit coupling) were worked out 
by Lax and Hopfield and by Streitwolf with single group 
theory . 98 ! 99 In this section we present a direct and de- 
tailed application of double group theory to derive spin- 
conserving and spin-flip matrix elements along arbitrary 
spin orientation directions. 

Rendering double group theory to explain the /- 
process spin flips is particularly simple since electron 
states involve only the two-dimensional (2D) irreducible 
representation (IR) Agf^i and no perturbation between 
other symmetry states is needed. For convenience, the 
/-process is studied by changing the spin orientation (n) 
while fixing the scattering to be between the +x and +y 
valleys. While the total number of independent (scat- 
tering) constants is not changed, the matrix elements 
depend on the chosen spin orientation. As such, the 
character table is not sufficient in mapping the depen- 
dence on the spin orientation. We bypass this limitation 
and use the explicit IR matrices. In this approach spin- 
conserving (in momentum scattering) and spin-flip tran- 
sitions are treated on equal footing. In the rest of this 
section we first study the case of spin orientation along a 
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FIG. 3: (Color online) (a) Phonon dispersion curves in silicon. 
Phonon modes at the vicinity of the red circles (black squares) 
take part in spin (momentum) mtervalley scattering. The /- 
process involves the zone-edge phonons near the E-symmetry 
axis (connecting the V and K points). The g- process involves 
phonons near the A-symmetry axis (between F and X points) . 
It will be shown that leading contributions are associated with 
phonon modes with the solid symbols. Intravalley scattering 
relates to all branches of long- wavelength phonons around the 
r point, (b) Relating phonons with Aj and E; symmetries to 
representative intervalley electron transitions. 



TABLE III: Non-trivial relevant IR matrices in a /-process 
between +x and +y valleys [E evaluated at ( — ko,ko,0), A 
at (fco, 0,0) or (0,fc ,0), k = 0.85 x 2rr/o]. Xa\ 2 = (x%T ■ 
Also shown is the effect of exchange operations on A star. 
Basis states in D1/2 is along ±2 in spin space. The || and _L 
directions in the mode subscript are relative to z. Eq is in 
unit of meV and Ds iS of meV • 27r/a. 
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valley axis. Then we will remove this restriction and pro- 
vide general expressions for the /-process spin-flip matrix 
elements [Eqs. (jlTa) - (fT7ij) and Table |P7j. 



A. Spin orientation parallel to a valley axis 

The /-process between opposite spin states in the +x 
and +y valleys depends on whether n is perpendicular 
to both valleys (n||z) or parallel to one of them (n||x or 
n|jy). For either case, we should consider four phonon 
symmetries represented by ID IRs, £1-4. These sym- 
metries dictate how each of the phonon-induced interac- 
tions, [H cp in Eq. ([5]) with Ej phonon], transforms 
under symmetry operations of the E group. In the left 
part of Fig. [3ja) we designate the phonon symmetries 
with the six phonon branches and in row 4 of Table IIIII 
we associate these symmetries with phonon modes. Row 
5 of the table lists the energies of these modes at the 
phonon wavevector that connects the valley centers. 

The nonvanishing matrix elements (k^l-HsJki) 
are obtained from selection rules that involve 
very few transformations ! 98 ' 100 These transfor- 
mations include common symmetry operations, 
g c G {(e|0), (e|0), (p z \r), {p z \r)}, of the little groups 
at ki, — k2 and q. They also include operations 
9e G {{pxy\0), (Px V \0), (faxylr), (^2s»|t)} that exchange 
ki and k 2 . The bar over operations denotes an addi- 
tional 2ir rotation. The character table of the nontrivial 
operations is shown at the upper three rows of Table HTT1 
The number of independent real constants involved 
in the matrix elements for each phonon symmetry is 
provided by consideration of these operations and time 



reversal symmetry, 



1 

2h 



9c 



{9c)x%{9c)xU9c) 



J2^M)xl(9e) 



(10) 



where ho = 4 is the number of g c or g e operations and 
XA 6 = XAi x Xi/2- The second sum in Eq. (fT0|) repre- 
sents the effect of time reversal symmetry and the minus 
sign takes into account the parity from the spinor basis 
and interaction H cp (sec, for example, Refs. jl0l| and 
}l02| for more details). By straightforwardly plugging 
the characters of Table IIIII into Eq. (fT0|) we get 



A/" El = 2, A/s 2 = 1, A/s 3 = 1, A/s 4 = 0. 



(11) 



Next, the interaction matrix elements 
(k 2 , s 2 |i?s 4 |ki, Si) between specific spin species are 
expressed in terms of Afx t independent scattering 
constants. To reach this goal, we first write down the IR 
matrix Da 6 for a given spin orientation. We choose the 
spin orientation conveniently along z direction, and then 
Da 6 = Da! x -D1/2 using Table ITTlU ^ Via appropriate 
group operations (see details in Appendix [A"t , we obtain 

(ka.ft, Iffsjkx,^) = D Sim +iDz lS , (12a) 
Qc 2 Az l^lkx,^) = D E2S - iD^ s , (12b) 
(k 2 ,^ \H Ss \ki,f z ) = Dx 3S +iDx 3S , (12c) 



and 



(k 2 ,4n \H St \k u H r 



(k 2 ,4n |flE 4 |kl 4n)*, (13a) 

-(k 2 ,^„ \H^\ki,)U)*,(13b) 
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TABLE IV: \Mj: i (ki, -ft" n ; k2, ii-a)/Ds iS \ for /-process spin flips between +x and +y valleys. The values of the scattering 
constants Ds i3 are given in Table [Till For each of the non- vanishing modes, Ei, the relative amplitude is provided for spin 
orientation (n) along any of the inequivalent high-symmetry crystal directions. Results between other valleys can all be obtained 
by trivial symmetry arguments. 
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up to a phase freedom for each matrix element. The scat- 
tering constants are all real numbers with D^ im much 
larger (about three orders of magnitude) than the rest. 
m denotes momentum and s denotes spin, with the rea- 
son more obvious in the general spin orientation case. 
Eq. ([13")) holds for all phonon modes and spin orienta- 
tions, and it expresses the effects of time reversal and 
space inversion (diamond structure). 

When the spin orientation is parallel to the axis of 
one of the valleys that participate in the /-process, then 
instead of changing the Aq basis states and the matrix 
form, we write the new spin states in terms of the pre- 
vious Ag basis (we call it the 'original basis'). As an 
example, for spin orientation along the x valley we use 
the rotation matrix [1, —1; 1, l]/\/2 as a unitary transfor- 
mation matrix and get^Si 

: (|k,^) + |k,^,»/V2, 

: (-|k,t,) + M,»/V2. (14) 

Using Eqs. (fT2j) and (|13|) we can obtain the results for 
spin orientation along one of the involved valleys (n || x; 
the case of n || y is equivalent). For Si, 



|k,t, 

Ma 




D 



S i rri 




(15a) 

Repeating the same procedure for the spin-flip case and 
for other phonon modes, we get 



B. General dependence on spin orientation 

So far we have restricted the spin orientation during 
an /-process spin flip to be along a valley axis (main 
crystallographic axis) . Removing this restriction adds to 
the anisotropy of spin relaxation processes and it allows 
one to make a direct comparison to a wide range of spin 
injection experiments. 

We generalized the previous derivation to an arbitrary 
spin orientation direction. We define n in terms of polar 
and azimuthal angles 9 and <f> with respect to the +z and 
+x directions. The new spin states relate to the original 
ones by an 'active' rotation matrix in spin space^i^ 

|k,fl„) ~ cos^|k,ts>+sin^|k,iy, 



■sin-e-^|k,^)+cos-|k,^). (16) 



Using the explicit matrix for interaction between 
|k, ft (-IJ-)z) basis states, as the example in Eq. (|15ap . the 
spin-flip matrix elements under the new spin orientation 
are 



(k 2 ,frn|ffE 1 |k l! fr n ) 

<k 2 ,^n|i?E 1 |k 1 ,fr n ) 

(k 2 ,frn|^s 2 |k 1 ,fr n ) 

(k 3j ^n|fl'EJkl I ftn) 



(k 2 ,4x 


\H Sl \k u 1t x ) = -iD Sls , 


(15b) 




(k 2 ,fr* 




(15c) 


(k 2 ,fr n |ff£ 3 |ki,frn 


(k 2 ,4x 


\H S2 \k u i[ x ) = D^ 2S , 


(15d) 




(k 2 ,fr* 


|i?s 3 |ki,1rx) = iDz 3S , 


(15e) 


(k 2 ,4n|ffs 3 |ki,fr„ 


(k 2 ,4x 


Iff^lkx,^) = % s . 


(15f) 




Together with Eq. 


[13)). the set of expressions in Eqs. ([12]) 




and (|15|) complete 


the results of inequivalent types 


of/- 





Dv im + i cos 6»D Sl 
-isin0e^£> SlS , 

smdsm((b + 



(1 + i) sm —e 



+ (l-i)cos 2 - 



(17a) 
(17b) 

Dz 2S , (17c) 



(17d) 



37T 

iv / 2sin6»sin((/) + — )£>s 3S , (17e) 
(l-i)sin 2 -e 2 ^ 



+ (1 + i)cos 2 



Dv 3S , (17f) 



process scattering when the spin orientation is set along a 
valley axis. Values of the scattering constants, Ds iS , are 
obtained by numerical calculations and they are listed 
in the last row of Table IIII1 Preferably, their values are 
extracted from experiments which cover some parameter 
range (e.g., temperature, stress and external fields). 



with Eq. (fl~3"|) for other matrix elements. Table IIVI lists 
the relative magnitudes of the squared spin- flip matrix el- 
ements for n along all inequivalent high-symmetry direc- 
tions of the crystal (they are often the spin orientations of 
injected electrons). Analysis of the spin relaxation time 
due to /-process spin flips will be given in Sec. IVIII 
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V. INTRAVALLEY AND INTERVALLEY 
^-PROCESS SPIN FLIPS 

In this section we will present a rigorous procedure to 
reach at compact intravalley and g-process spin-flip ma- 
trix elements. Before embarking on the theory, we discuss 
key considerations that underly the analysis and exem- 
plify their outcomes via representative results (Table |V} . 
This choice allows one to understand the most impor- 
tant physical parameters without delving into details of 
the analysis (which are provided in Sees. IV Al through 
EE). 

Spin-flip matrix elements vanish for intravalley (g- 
process) scattering if wavevectors of the initial and final 
states are the same (opposite). The wavevector power- 
law dependence was mentioned in Table U along with 
the fact that suppression of zeroth-order terms in these 
processes leads to relative slow spin relaxation compared 
with the /-process. However, intravalley and g-process 
become more important when shear strain is applied 
since different valley minima are split thus suppressing 
the scattering by the /-process i 15 ' 22 The intravalley spin 
flips are also important at low temperature due to the 
larger population of long wavelength acoustic phonons 
[smaller energies; See Fig. (3ja)]. 

It is instrumental to compare spin flips with momen- 
tum scattering of which the study is more established. 
Following the ingenious connection with deformation po- 
tential parameters by Bardeen and Shockleyji^li Herring 
and Vogt derived a detailed angle dependence of intraval- 
ley momentum scattering due to interaction with acous- 
tic phonon.-^£ In spin flips, the connection with deforma- 
tion potential is more subtle and complicated by the de- 
pendence on high-order wavevector components. Later, 
we will derive explicit forms while making no a prior as- 
sumptions about the form of possible deformation poten- 
tial parameters. This approach allows us to identify the 
crucial role of the coupling between the lowest pair of con- 
duction bands in setting the intravalley spin relaxation. 
The dependence of spin-flip matrix elements on high- 
order wavevector components makes this coupling effec- 
tive in spin relaxation (while being marginal in momen- 
tum relaxation). The relaxation rate becomes inversely 
proportional to the square of the energy gap between the 
conduction bands at the valley center, Ac w 0.5 eV. 

Implications of the above considerations are mani- 
fested in the spin-flip expressions of intravalley and g- 
process scattering. The matrix element of intravalley 
spin flips will be shown to consist of four factors that 
represent different aspects of the above considerations 



|M^ tra (k! ^ n ,k 2 4n)| 2 



A c 



DlS n (q)h(q) 



The factor |7y|/Ac « 0.03 A is calculated from band struc- 
ture parameters that originate from spin-dependent kp 
perturbation terms in the Hamiltonian. The second fac- 
tor, D\, depends on the phonon mode. Scattering with 



TABLE V: Squared spin-flip matrix element Af s f j 2 , induced 
by all types of phonon modes in intravalley scattering, and by 
the LA phonon mode in g-process scattering, q = k2 — ki and 
K = (k2 + ki)/2. Valley centers are set along z axis. Spin 
orientations n are taken along all inequivalent crystal sym- 
metry directions. Results in other valleys can be obtained by 
trivial symmetry arguments. See text for related parameters 
and further explanations. 
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long-wavelength acoustic modes is governed by a defor- 
mation potential constant, Dta/la = D' xy , that couples 
the two lowest conduction bands. Its explicit integral 
expression will be given in Eq. (|46[) and its numerical 
solution (via EPM calculation) yields a value of 6 eV. 
This value agrees well with the measured energy splitting, 
4D' xy e xy , of the lowest conduction bands at the X point 
when applying a shear strain i ° 8 ' 109 Scattering with long- 
wavelength optical modes is governed by the constant 
Dto/lo = D op rj 5 eV • 2ir/a which originates from the 
out-of-phase motion of atoms in the primitive cell. This 
parameter is also calculated from the coupling between 
the lowest pair of conduction bands [Eq. (JUJ)]. The re- 
maining two factors in the above expression depend on 
the valley position. Their forms for electrons that reside 
in the z valley are listed in Table [V] I\ and S n reflects 
the dependence on phonon properties and on the electron 
spin orientation, respectively. The right column of the ta- 
ble lists the values of I\ for all long-wavelength phonon 
modes and the left column lists the values of S n for five 
inequivalent high-symmetry directions in the z valley. 

The g-process spin-flip will be shown to share several 
properties with the intravalley case. For electron transi- 
tion between the ±z valleys, the dominant g-process spin- 
flip mechanism will be shown to originate from scattering 
with LA phonon modes, 

|Mf A (k! ^ n ,k 2 ^ n )| 2 =£ 2 s S n (K), 



where its prefactor D„ 



O.leV is dominated by the 
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sum of two scattering constants. The first is related to 
the deformation potential of the lowest conduction band 
[Eqs. (|E12j) - (|E14j) ] and the second couples the lowest con- 
duction band with upper valence band by spin-orbit mod- 
ulated clcctron-phonon interaction [Eq. (IF12I) ]. 

The rest of this section is organized in the follow- 
ing logical order. In Sec. IV Al we derive selection rules 
that pertain to both the band structure Hamiltonian 
and the electron-phonon interaction. The obtained elec- 
tron state vectors are used in Sees. IV Bl and IV CI where 
we present the core derivation for intravalley and g- 
process spin-flip matrix elements. In these parts we de- 
rive exact forms of the various scattering constants. Fi- 
nally, the dependence on spin orientation (n) is given 
in Sec. IV Dl Readers who are not interested in the full 
derivation may directly skip to the most general expres- 
sions in Eqs. (|45p - (|48|) and Fig. [5] for intravalley spin- 
flips in the z valley with n |j z; in Eqs. (|54"l) -([55 j) and 
Fig.[5]for g-process spin-flips between the ±z valleys with 
n || z; and in Eqs. (|57)) - ([5^)) for general spin orientations 
(fi = [cos 4> sin 9, sin (f> sin 9, cos 9] ) . 



A. X-point selection rules and spin-dependent 
eigenstates of the k ■ p Hamiltonian 

In this part we present group theory results and an- 
alytically quantify the signature of spin-orbit coupling 
on electronic states. Since the bottom of the conduction 
band is at the vicinity of the A point [see Fig. 0Ja)], 
we employ a compact k • p model with a small set of 
basis states that pertain to the symmetry of this point. 
Findings of this model will be used to derive spin-flip ma- 
trix elements. These findings will also be benchmarked 
against numerical calculations of an EPM that includes 
spin-orbit coupling.— 

For intravalley scattering, working with A-point space 
group is as effective as working with the group of the A 
axis (position of the valley center), and bares favorable 
features over the latter choice-^ The initial and final 
states of a ^-process scattering can be expanded to the 
same X point, making it possible to relate the matrix 
element with a deformation potential parameter, though 
with modifications (to be discussed). 

We first mention spin independent featuresi^S 8 . Due 
to the symmetry of the crystal, only one of the six 
conduction band valleys is studied and we arbitrary 
identify it as the valley along the +z crystallographic 
axis. We choose a basis of four eigenstates of the 
Hamiltonian Hq = p 2 /2rriQ + V{r) and denote them by 
{\X$),\X{),\X?),\X%)}, where Vx(r) = (r|A). These 
A-point states are associated with the lowest pair of 
conduction bands and upper pair of valence bands at 
k = (0, 0, kx ), where kx = 2ir/a [see Fig.@Ja)]. Inclusion 
of the A4 valence states is imperative since they bring the 
mass anisotropy of conduction electrons. The Ai and 
A4 nomenclature denotes 2D IRs of the space group G| 2 
which describes the symmetry of the A pointi 108 ' 111 " -^ 




r X U u (2rc/a) 0.02 



FIG. 4: (Color online) (a) Energy band structure along the 
A-axis between the X = (0, 0, l)2n/a and T — (0, 0, 0) points 
and along the E-axis between X and U = (1/4, 1/4, l)2-7r/a 
points [see Fig. [TJa)] . Solid lines are the results of EPM and 
dash lines are taken from Eq. (|C4[) of the Appendix. Ex x — 
Ex 4 = E g ,x and we set Ex x = 0. Ac is the energy gap 
between the two conduction bands at the valley center ko — 
(0, 0, 0.85)27r/a. (b) Conduction band dispersion along the 
^-symmetry axis between X and W = (0, 1/2, l)27r/a points. 
Along this direction the two-band conduction degeneracy can 
only be lifted by spin-orbit coupling (spin hot-spot), (c) The 
induced energy splitting by zooming-in closer to the X point. 



The dimensionality complies with the two-band degen- 
eracies at the A point of diamond crystal structures 
(due to time reversal and glide reflection symmetries). 
The indexing of the basis states implies of their com- 
patibility relations when going from the A point into 
the A axis ({|A y ), |Ai), |Af>, |A|)})4fia For example, 
ipA 2 ,,k z { r ) — e % ( k *~ kx ^ z ip X 2< (r) . Finally, we mention 
that for the following derivations, the symmetry prop- 
erties rather than explicit functional forms of the basis 
states are important. 

When adding the spin degree of freedom the new ba- 
sis set reads { \X?,f), \X?,±), \X^), \XH), |Af ,t), 
|Af,4.), [A']','!"), |A|,4,) }. Two remarks on this basis are 
in place. First, inclusion of the A4 states in the basis set 
enables us to capture the salient spin-dependent prop- 
erties of conduction electrons (originate from spin-orbit 
coupling with the upper pair of valence bands).— Sec- 
ond, the spin-orbit coupling is treated as a perturbation 
and, accordingly, the new basis states are not eigenstates 
of the spin-dependent Hamiltonian at the X point. In 
Eqs. ([26 |) -(f27 | of this section we find specific forms of 
these states for electrons around the bottom of the +z 
conduction valley for spin orientation n||z [see Eq. flSJ)]. 
In sections IV Cl we discuss the needed changes when con- 
sidering electron states from opposite valleys. In Sec- 
tion [VD] we show the anisotropy of spin- flip processes by 
varying n. 

Group theory is invoked to construct the band struc- 
ture Hamiltonian matrix and matrix elements of electron- 
phonon interactions . 101 i 114 i 116 i 117 Both depend on cou- 
plings of Xi basis states, of A 4 basis states, and of A x 
with X4 basis states. Using the character table of G§ 2 
(Appendix |BJ), these couplings are clarified by the follow- 
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ing decomposition: 



,118 





Xi a 


3X X 




X 4 a 






X 4 <i 


5Xi 



Mi © M4 © M' 2 © M's , (18a) 
M l ®M 4 ,®M[®M' 4i , (18b) 
M 5 © M' 5 . (18c) 

M5 and M5 are 2D IRs and the remaining Mj arc ID 
IRs. Next we relate these decomposition with specific lin- 
ear combination of state productsJ^i For example, one 

of the X\ states coupling [Eq. (|18aj) ]. ip* 2 , (r)tp x i (r) + 

x x 1 

ip X 2> (r)ip* x i (r) is found to transform as M 4 (apply the 

operations of G 2 2 on it). We associate this linear combi- 
nation with [0, 1; 1, 0], or p x , in H cc block (In this paper, 
Pi and cr, identify, respectively, Pauli matrices in the 2D 
IRs product and spin space). We find the following asso- 
ciations for product combinations in intraband coupling, 

H cc : Mi-H-X, M A ^p x , M' 2 ^p y , M'^p g , (19a) 
H vv : Mi+*X, M 4 ^p x , M[^p y , M'^p z . (19b) 

In case of interband coupling where the decomposed IRs 
(M5 and Mg) are 2D, we classify pairs of product combi- 
nations. For example, {ipx^x 2 ' + ^x^^xl 1 0xf V'xJ + 
^ s ^2'} belongs to M 5 . We find 



M 5 ^{Py,p Z }, M' s {1 , p x } , 



(19c) 



Non- vanishing matrix elements of any interaction (i.e., 
selection rules) are readily identified once we specify how 
terms in that interaction transform under symmetry op- 
erations. These terms usually can be made into parts 
that transform as components of vectors (r) or axial vec- 
tors (R). The longitudinal component of the vector (ax- 
ial vector) transforms as M' 3 (Ms), and transverse com- 
ponents as M 5 (M5): 

z ~ M3, {x,y} ~ Af B , R z ~ M 3 , {R x ,R y } ~ M 5 .(20a) 

where the longitudinal component is along the valley axis 
(z in our choice). Besides, the symmetric and antisym- 
metric potentials [Eq. ([7])] are even and odd under space 
inversion and transform as 



V 4 



Mi 



V- 



(20b) 



By the fundamental theorem of group theory, a non- 
vanishing matrix element results when the interaction 
operator and the states product decomposition belong to 
the sameii^ IR(s). The integrand of the non- vanishing 
matrix element as a whole then belongs to the identity IR 
(in this case, Mi). From Eqs. (|18c[) and (|20aj) . the trans- 
verse components of a vector (M5) or axial vector (M5) 
can couple conduction band with valence band. How- 
ever, since Mi only appears once in either M5 © M5 or 
M' b ®M' bl only one of the linear combinations of the oper- 
ator and state product can survive and the unique linear 
combination read 



for M 5 ©M 5 or M^©M^, respectively. The IR product de- 
compositions and the IR assignments with (axial) vector 
are derived from G\ 2 character table (see Appendix IB"]) . 

Transformation under time reversal operation (7~) is 
equally essential in obtaining selection rules (in addition 
to spatial operations). Time reversal symmetry connects 
matrix elements by^Si 



{Xj\0\Xi) = (TXilTO^T-^TXj 



(22a) 



Physical operators O have different parities under time 
reversal operation, 

r^r, W^W, p^-p, Wxp -> -Wxp,(22b) 

and this feature distinguishes between vectors such as 
Vy and p. T acting on \Xi), by our convention, ex- 
changes the basis states?^ 



IX^oIXi 1 ), |Xf)o|Xf> 



(22c) 



Therefore, for the diagonal blocks (H cc and H vv ) 

X^t-X, p x -> p x , p y -> p y , p z -> -p z (22d) 



under time reversal operation in light of Eq. (|22a[) . A 
matrix clement vanishes if under time reversal operation 
the changes of sign are opposite for the state product 
[Eq. (|22dj) ] and the interaction [Eq. ([22b)) ]. 

The spin-dependent k ■ p Hamiltonian matrix can be 
prescribed using the above group theory analysis, with 
perturbation Hamiltonian (H = Hq + Hi) 

Hl = (¥JfL + Hk ' ■ P) X + W) x (P + ^)l' ff )(23 ) 



2mo 



m 



4toqC 2 



where k' = (k x ,k y ,k z 



- kx ) is measured from the X point 



and V(r) is the crystal potential Ylj V+( r ~ R-j) 
Hamiltonian matrix blocks are thus, respectively 



H\cc — 

Hlvv — 
Hive = 



2mo 

h 2 



(k' 2 l®l + 2k[ ) k' zPz ®l), 
k' 2 I © I, 



2mo 

-iP(k x p y + ik y p z ) ©I + iA x (p x © o y 
-X®a x ) + a[k' z (ip z ®o x -p y ® a y ) 
+(kyp v - ik x p z ) © a 



. The 

(24a) 
(24b) 



(24c) 



ixpy - yp z — Mi and R X I - R y p x ~ Mi , 



(21) 



We stress that dV(r)/dz in H\ cc block vanishes con- 
sidering time reversal symmetry.The spin-independent 
constants in Eq. (|24[) are empirically known from 
experiments^^ They relate to the position of the con- 
duction band minima k = (0,0, k' + kx), and to the 
momentum matrix element that sets the mass anisotropy 
(P). The spin-dependent parameters (Ax and a) are 
calculated by the empirical pscudopotcntial model. Ta- 
ble IVII lists these and other parameters that we use in 
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TABLE VI: Parameters of bulk silicon. The lattice constant 
is a=5.43 A, the X-point energy gap is E g , x ~ 4.3 eV, and 
the free electron mass is mo. 



k' 


-0.15 


2-k/o, 


Eq. (|Cla|| 
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eV-a/2vr 


Eq. ||Clb| 
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meV 


Eq. \C\c\ 
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eV 


2h 2 k' 2 /m 
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16.7 
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r) = 2iPAx/E g , x 
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18.9 
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rf = 2iPA' x /E g , x 



this paper. In Appendix [Cj we provide their integral ex- 
pressions and we also apply a partitioning technique to 
analytically diagonalize the total Hamiltonian matrix 



H\cc H lvc 

Hive H\dd E, 



y.x 



(25) 



where E g x is the energy-gap between the conduction and 
valence band at the X point [see Fig. Ufa)]. 

Appendix [Dl provides a general procedure to align the 
resulting degenerate states to spin states that satisfy 
Eq. ((6J. For spin orientation along the valley axis (n||z), 
degenerate states in a conduction band valley read 



|k,4) 



exp(ik' • r) [A(k)| f) 
exp(ik' t) [A(k)||) 



B(k)|i)]|X), 
B*(k)|t)]|X) (26) 



where 



|X) = 
A(k) = 

B(k) = 



\X*),\X}),\Xf),\X, 



2P 2 k r k 

■n 



1 



Eg,x£±C 

A~ c 



1- 



2E 2 

k z — k 



Pk x 
E~ 



E a 



ir)'(k x +ik y ) 
2E n 



E, 



Ac E, 



g,x 



'g,x 

E~ 



y.x 



Ac E, 



C ^g,x 



.(27) 



As shown in Appendix. [Ej momentum scattering is gov- 
erned only by the zeroth- and first-order wavevector com- 
ponents of A(k). For spin relaxation, on the other hand, 
the 'negligible' spin-orbit coupling coefficients become 
crucial [B(k) vector]. In writing the coefficients of A(k) 
and B(k), we have kept only those terms that are rele- 
vant to the spin relaxation analysis of the following sec- 
tions. Terms that scale with Ac / (2E gX ) « 1/17 are 
omitted due to their negligible effect. The energy gap 
between the two conduction bands at the valley center 
Ac ~ 0.5 eV [see Fig.UJa)] will be extensively used when 
dealing with intravalley spin relaxation. The main differ- 
ence of the above solution from Ref. [2(| is that we have 
identified all quadratic-wavevector terms that play im- 
portant roles in intravalley spin-flip processes. We have 



checked the accuracy of these expressions with numeri- 
cal EPM solutions The main difference between them 
comes from the omission of the lower valence band states 
[lower A"i point in Fig. Eta)]. These states are respon- 
sible for making the longitudinal mass slightly less than 
too- In spite of missing some k z dependent components, 
the error amounts to a few percent and is irrelevant in 
intravalley spin relaxation (to be shown later). 



B. Intravalley spin flips 

In this subsection we present the theoretical procedure 
for deriving intravalley spin-flip matrix elements for all 
types of phonon modes. Scattering with long-wavelcngth 
TA phonon modes will be shown to dominate the in- 
travalley spin relaxation at low temperature. This prop- 
erty is largely set by the phonon dispersion around the 
r point. Figure 02a) shows that the energy of long- 
wavelength acoustic phonons is linear in q with a smaller 
slope for TA, and approaches zero at T point. For long- 
wavelength optical modes, the phonon energy is almost 
wavevector independent («63 meV). However, when the 
temperature increases, optical modes become more im- 
portant and together with the LA mode they provide 
about half of the intravalley spin relaxation at room tem- 
perature. The reason is twofold. Spin-flip matrix ele- 
ments due to scattering with optical modes are linear 
in q while being quadratic with acoustic modes (to be 
proven below). Second, the rise of phonon population 
with temperature is much faster for optical modes. 

In Appendix [El by invoking selection rules we have 
redcrived the leading-order matrix elements of intraval- 
ley momentum scattering, which is typically described 
in terms of deformation potential constants ! 106 ' 107 When 
extended to analyze the intravalley spin-flip matrix ele- 
ments, this group theory approach will be appreciated by 
its ability to reach essential terms efficiently. It also pro- 
vides flexibility in the sense that no a priori knowledge of 
the connection between deformation potential quantities 
and spin-flip processes is needed. 



long-wavelength acoustic phonon modes 

Spin-flip matrix elements due to electron interaction 
with long-wavelength acoustic phonon modes have a 
quadratic-wavevector dependence. We show it by start- 
ing with the most general form, 



xl k 1= K-a f ; k 2 =K + SjL) = (28) 



M A (k 1= K-a ft ; k 2 =K 
(K+f A |tf(q) ■ VV+(r) +£-(q) ■ W_(r)|K-f ,ft) 



Mgf\ denotes terms of j th order in q due to A phonon 
mode. In the next step we expand the states in Eq. (f2"8"]l 
in increasing orders of q. We show how the zeroth- and 
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first-order terms vanish {M^\ = M^' x = 0), and then we where for the case in hand space inversion provides 



(2) 

find the form of the dominant quadratic terms (M f ) . In 
this wavevector-order analysis, the q dependence of the 
in-phase and out-of-phase polarization vectors has to be 
taken into account. In Table ITTI we showed these depen- 
dencies in the long-wavelength regime which is relevant 
for intravalley scattering (q <C 2ir/a). The out-of-phase 
vector, £j"(q), * s l mear m Q while £iT(q) has a zeroth- 
order dependence (e.g., qi/q terms). Hereafter, we ab- 
breviate the notation £*(q) and as ^ ± and M^f . 

Zeroth-order: The lowest-order spin-flip matrix ele- 
ment has the form 



M 



(0) 



(K,^|£+- W+|K,ft). 



(29) 



We exclude the out-of-phase part since the acoustic inter- 
action at the infinite wavelength limit is governed solely 
by the in-phase motion of atoms in the primitive cell 
[£~(q = 0) = and |£+(q = 0)| = 1]. In this limit, the 
phonon-induced interaction reduces to displacement of 



the entire crystal, J2j ^rV+(r 
on the relation 



Wcr 



crystal 



i[p,H]/H, 



VV cr v S tai- Based 



(30) 



the coupling of the in-phase part between spin-degenerate 
eigenstates of H vanishes. When breaking this vanish- 
ing matrix element into parts that come from interaction 
with the bare symmetrical potential and with its spin- 
orbit coupling part, 

(KA |W+|K,f) + (KA |Wf°|K,ft) = , (31) 

the sum is zero but each of these two contributions is 
finite. This result was first pointed out by Elliott^ Later 
we will make use of this property 

First-order: To write the spin-flip matrix elements of 
this order we begin by linearizing the state, 



|K±f )S ) 



- |K, S )±iq.£|K, S )+0(<7 2 ), (32) 
where C denotes the derivative of the state in k-space, 
\h + Ski,s) - |k,s) 



£Jk, s) = lim 

<5k^0 



8k l 



(33) 



Substituting Eq. (|32[) into Eq. (f28j) and considering the 
wavevector dependence of ^ ± , the first-order spin-flip ma- 
trix elements are, 



M 



(i) _ 



si' 



= (KA\C -W_|K,ft) 



(34) 



+ i(K,4| (q-£ f ) (£+■ W+)|K,f) 

- i<K,^|(£+-W + )(q-£)|K,ft) . 

Combining space inversion and time-reversal symmetries, 
the out-of-phase term can be shown to vanish (first line 
on the right-hand side) , and the in-phase terms to cancel 
each other (second and third lines). These facts follow 
the relations 



S\kA) = |-k,ft) 
and time reversal provides 



5|k,^) = |-k,4) 
SCS' 1 = C , (35) 



T|k,ft> = |-k^) , r|k,^) = -|-k,f) 
rwlr- 1 = w± , r{c^j- x = c . (36) 

Yafet separated the in-phase terms into two parts to 
prove them vanishing J^- Together with Eq. (|3ip . it is the 
celebrated Elliott- Yafet cancellation. We emphasize that 
the first-order (lincar-in-q) Elliott and Yafet terms van- 
ish separately, rather than interfere destructively. All 
in all, the zero- and first-order spin-flip matrix element 
identically vanish, M s ( f 0) = Afjp = 0. 

Second-order: It is the lowest order at which spin-flip 
matrix elements due to electron interaction with acoustic 
phonon modes do not vanish. At this order, states are 
expanded by 



|K± 



|K, s )±iq-£|K, s ) 



±q® 2 -£® 2 |K,s), (37) 



where the vector components of C were formally defined 
in Eq. ([33]), and q® 2 ■ C® 2 denotes the scalar product of 
two second-rank tensors. An explicit form of this state 
was derived by a spin-dependent k • p expansion of the 
X-point basis states [Eqs. (|2l)|) - (|2"7) ]. Using this basis, 
the general spin-flip matrix element is converted to 

M (k x = K - a t, k 2 = K + f,4) = Alf + 0(q 3 ) w 



(X /1 |e- i * I [r+ I/ (k 1 ,k 2 )+r- !/ (k 1) k 2 )]e-^ £ |X I/ ). (38) 



r^ : 1 ,(ki,k2) arc scalars formed by products of column 
and row spinors with the potential matrix, 



r± (ki,k 2 ) 



(39) 



[-S XM (k x ),A XM (ki)] (£± ■ W±, K (r)) 



^„(k 2 ) 



where the coefficients of A(k) and B(k) are taken from 
Eq. (|2"T)) , and the K-dependent potential has the form 



V±,K = 



V±l- 
V±T- 



h 

Arrive 2 



[VV± x (p + KK')} ■ a . (40) 



The bare potential, V± , is diagonal and generates Elliott 
products of the type A^B V . The spin-orbit coupling po- 
tential, V±° k , generates Yafet products of all types but 
the dominant signature comes from A m A n terms due to 
the smallness of the coefficients in B(k). 

We identify the general forms of second-order matrix 
elements between the X-point basis states. After ex- 
panding the exponential and T terms in Eqs. (|38l) - ([3"9")) 
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into power series, quadratic terms in q are classified by 
six integrals that read 

M S ] = EEw, (4i) 

fi,v n—1 



;5 



2 \ X V 



= x 



X, 



x, 



k=k 



k=k 



C M (k){qT,r-VV_, k (r)}C„(k) 
q • V k C M (k)) (C ■ VV_, k (r)) C„(k) 

q • V k C M (k)) (€+ • VV+, k (r)) (q • V k C„(k)) 
q ■ V k C, t (k)) {q ■ r, £+ ■ VV+, k (r)} C„(k 
q^ 2 • Vf C M (k)) (|+ • VV+, k (r)) C„(k 
C M (k)q-r(e+- VV + , k (r))q.rC„(k) 



A', 



k=k 



X„ 



k=k 



k=k 



A', 



k=k 



(42a) 
(42b) 
(42c) 
(42d) 
(42e) 
(42f) 



r 



where {£/, W} = UW + WU and 

C v (k)= [A Xu (k),B Xu (k)] T , 
C M (k) = [-B x ^k),A X)i (k)] . 



(43) 



In writing these matrix elements we have used time re- 
versal and space inversion to unify terms. The mean 
wavevector K is replaced with ko, which brings approxi- 
mation of 0(q 3 ). 

At this point, the integrals of Eq. (|4"2"j) can be numeri- 
cally calculated in a straightforward manner. Given the 
four basis states and the above integral forms, this nu- 
merical procedure involves a few hundreds of different 
types of space integrations. Performing selection rules for 
all of them also demands a large amount of work. How- 
ever, we can avoid these labors and come to an accurate 
and compact matrix element [Eq. (|45[t] . This possibility 
is enabled by three observations that greatly simplify the 
analysis of Eq. (|42l) . These observations will also allow 
us to connect the derived matrix element with an exper- 
imentally known deformation potential quantity. 

Observation 1 : Integration of bare potential re- 
lated terms such as (X fJi \W±\X„) , (X^n W±\X U ) 
and {Xf^TiTj W±\X V ) yield numbers of the same or- 
der of magnitude given that fj is measured in units 
of a/2n. Moreover, these numbers are comparable for 
both the in- phase (+) and out-of- phase (— ) parts. The 
same applies to the integrals of the spin-orbit cou- 
pling potential, {X^V^ M \X V ), (X fJi \{r l ,'VVl k }\X ly ) 
and (X f j l \ri'W^° k rj\X v ). These conjectures are backed 
by explicit numerical calculations. The physical rationale 
is that these potentials are significant within the size of 
a primitive cell where r, is of the order of 1. In what 
follows we keep consistency and measure length in units 



of a/2n. Accordingly, band-structure parameters such as 
P, r\ or a are approached as energy scales, and compared 
directly with E g , x or A c (see Table ED . 

Observation 2: The number of matrix elements is 
greatly reduced by estimating the amplitude of their 
coefficients. The largest Elliott and Yafct products 
scale, respectively, with \t]\/Ac (or \-q\P / AcE gX ) and 
2P 2 /(E g ^ x Ac)- For both cases, other products are 
significantly smaller. These products relate, respec- 
tively, to the bare and spin-orbit coupling poten- 
tials. These amplitudes are evaluated by inspection 
after substituting the explicit coefficients of A(k) and 
B(k) [Eq. dUD] into Eq. g2J. Notice that the terms 
{C(k), V k C(k), V 2 .C(k)}| k=ko are all constants without 
k dependence (and similar for C). 

Observation 3: Prior to the application of selection 
rules, contributions of Elliott and Yafet processes to in- 
travalley spin relaxation are conceivably comparable. To 
understand this physics we first recall the zeroth-order 
Elliott- Yafet cancelation [Eq. (|3ip]. For states at the 
valley center (k = ko) this cancellation means 



(XxIW-jAi) 



(AilW+IX^ 



A 



x 



(44) 



Taking into account observation 1 , we can generalize this 
order of magnitude to ratios between all Yafet-related 
integrals and all Elliott-related ones. Dominant Elliott 
processes relate to matrix elements of the type that ap- 
pears in the denominator of Eq. (|4"4"|) multiplied by r\j Aq 
(see observation 2). Similarly, dominant Yafet processes 
relate to matrix elements of the type that appear in the 
numerator multiplied by 2P 2 /(E gX Ac)- The overall ra- 
tio between Elliott and Yafet processes is thus of the 
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order of unity (~ P/E g x ). 

Application of observations 1-3 in Eq. (|42|) results in a 
handful of matrix elements that are worth examination. 
For example, the integral classes Ip,,u;i and I^^fi are elim- 
inated on ground of their small coefficients [A' x /E g x and 
1 rather than r\l Ac and IP 2 j[E g x Ac) in observation 2]. 
Therefore, 



and I, 



[i,v; 6 







for any possible basis state combination (X„ and X v ). 

We invoke group theory and evaluate the remaining 
integrals in Eq. (|42|) . Group theory results [Eqs. (fl"8|) - 
([22]) and their discussion] are extensively utilized, cou- 
pling integrals are expressed analytically by reasonable 
approximations, and the internal displacement of silicon 
structure are carefully accounted. The detailed proce- 
dure is shown in Appendix [F] We reach at the intravalley 
spin-flip matrix elements due to interaction with acoustic 
phonon modes, 

M^ik.A ; ^) = ^D' xy e xy ,M)i^-i<ly) 



ir,D' x 



■HI , 



A 



c 



(q x -iq y ) &£t„(q) + g v $L(q) 



(45) 



where A={TAi, TA2, LA}, with a physical deformation 
potential constant D' xy whose complete integral expres- 
sion reads 



-P 2 m 



■HI 



- x 



h 2 xyz 

A,B 



dz 

dV at (? - r a ) 



(46) 



d.v 



x\ 



6 eV. 
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FIG. 5: (Color online) Analytical (solid lines) and numerical 
results (-x- lines) of spin-flip matrix element as a function of 
q = ki — k2 along five different scattering directions. Each 
direction is represented by a different line color. The electron 
resides in the z-valley and its spin is orientated along the val- 
ley axis. The analytical results for TA, LA, TO and LO modes 
are taken from Table [V]for n || [001]. Numerical results rely 
on empirical pseudopotential and adiabatic bond charge mod- 
els (to be discussed in Sec. I VILA")) . Analytical and numerical 
curves agree very well for all cases where the quadratic depen- 
dence of the acoustic cases holds best for small q values. Note 
that analytical results for different q directions may overlap 
each other. 



where x={TOi, TO2, LO}, with the associated scatter- 
ing integral 



long-wavelength optical phonon modes 

Using the gained knowledge of spin-flip processes with 
acoustic phonon modes, we can readily derive the opti- 
cal case. The optical phonon modes have a dominant 
out-of-phase polarization vector, £~ p , with magnitude of 
about unity at the long- wavelength regime (Table ITTT) . 
From the detailed analysis of the out-of-phase acoustic 
phonon modes we could have recognized the fact that the 
antisymmetric interaction, VF_, has a leading-order ma- 
trix element that is linear in q. Separating the interaction 
with optical modes into in-phase and out-of-phase parts, 
we find in this case that the out-of-phase part dominates 
the in-phase by two orders of magnitude. This result is a 
consequence of the q dependence of the polarization vec- 
tors as well as of the interaction integrals. The leading 
integral term is simply 

M^ tra (k!^ ; k 2 ,4) = ^A>p (& - iq y )^,M , (47) 



D op = (X 2 ' \dV- /dz\X\) w 5 eV • 27r/a. (48) 

The ability to investigate for each individual mode is 
desired under conditions of anisotropic fields or stress. 
Then, the relative importance of different modes may 
vary depending on the symmetry breaking mechanism. 

We render the elastic continuum approximation in or- 
der to achieve complete and analytical q-dependent ma- 
trix elements and to facilitate analytical integrations. 
The polarization vectors, ^ ± (q), in Eqs. (|45|) and (|47|) 
are then replaced by the expressions of Table [TTJ We 
obtain |M sf | 2 for each of the modes (TA, LA, TO and 
LO)J^£ The final expressions are given in Table fVl (the 
case of n|| [001]). To examine the accuracy of the q de- 
pendence, Fig. [5] shows the analytical intravalley results 
of Table [V] next to numerically calculated matrix ele- 
ments along inequivalent high-symmetry scattering angle 
directions. The numerical procedure will be discussed in 
Sec. IVIII together with application of the derived matrix 
elements to (analytically) find the spin lifetime expres- 
sion. 
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C. Intervalley g-process spin flips 

g-process scattering angles are strongly directional 
(q = k 2 — ki ps — 2k ). The leading spin- flip matrix 
element of a g-proccss depends on the mean wavevector 
K = (ki +k2)/2 (and not on q as in the intravallcy case; 
See Table IJ). When k 2 = — k x the opposite spin states 
are Kramers conjugates and their coupling via scattering 
with any type of phonon vanishes*^ As will be shown 
below, the LA mode leads the <?-process spin relaxation 
where its matrix element is linear in K. The scattering 
constant, however, is not as concise as in the intraval- 
lcy case. This is conceivable because the translational 
factors, e ±zk °' r , of the initial and final states are to be 
expanded around a common X point (k' = 0). Other 
than the dominant LA phonon mode, we also discuss the 
general shape of the matrix element due to scattering 
with TA phonon modes. Their contribution to ^-process 
spin relaxation at room temperature will be shown to be 
non-negligible in comparison with the LA mode. This 
property is analogous to g-process momentum scatter- 
ing, where TA phonon modes with a higher-order ma- 
trix element but lower energy (compared to the lead- 
ing LO phonon) are important in describing the charge 
transport In what follows we continue to work with 
the X-point basis states which can be readily related to 
the A-axis basis states via i>Ai,k z — e lkzZ i^Xi- 

For ki and k 2 at the vicinity of the ±z valley center, 
we expand the matrix elements around K = 0, 

= ^(§+K,^|^ ± (q).VV ± |-§+K,t) 



,x 10 



r(l) 



0(K)0(\q/2 + ko|) + 0(X 3 ), (49) 

and explain each term separately. In a g-process £~ (q) 
is not treated as a small quantity compared to 1. The 
in-phase and out-of-phase interaction are treated on an 
equal footing. However, we will find that the out-of-phase 
contribution drops. 

Applying time reversal symmetry on a spinor state and 
k-derivative operator C [Eq. it is readily seen that 

the zeroth-order term (K = 0) vanishes 



M 



(0) 



± 



:f,4k ± (q)-VV±|-f,^)-0, (50) 



and that the two linear-in-K terms are equal to each 
other, 



(i) 



si 



= E<f^l(* ± - vv ±) ( K -£)l-§^> 

± 

£<f^ (^•K)(| ± -VV ± ) -lit). (51) 



When C acts on the translational part, e lk r , the rcsult- 
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FIG. 6: (Color online) Analytical (solid lines) and numerical 
results (-x- lines) of p-process spin-flip matrix elements as 
a function of K = |(ki + ka) due to scattering with LA 
phonon modes. Results are shown for five typical scattering 
directions where each is represented by a different line color. 
Analytical expressions are given in Table |V]and their results 
for K|| [001] and [100] overlap each other. Numerical results 
rely on empirical pseudopotential and adiabatic bond charge 
models (to be discussed in Sec. I VII A|) 



ing matrix element becomes 



M 



(i) 



sl'.l 



E<§ ^ I • VV ±' *K • r] | -§, ft) « 0, (52) 



where [U, W] = UW - WU . The Elliott part of this 
term naturally drops since VV and r commute. The 
Yafet part results in a small 'aK' factor and a weak cou- 
pling between basis states (i.e., non-dominant coefficient 

products). Thus, Mf, can be safely discarded. 

The second part of the first-order term, M^' 2 , results 
from operating C on the k-dcpcndent coefficient C(k) 
and C(k) [Eq. (gSJ] of the bra and ket states in Eq. ([ST]) . 
Time reversal symmetry is utilized again where following 
Eq. (J361) we get that (f,4J. | = (T| - f,fr)) f , and the 
action of T on individual X-point basis states follows 
Eq. (|22c|) . All in all, the second part of the first-order 
term reads 



2^^(l, K- VkC^koK*" 
(£± • VV±)a,(k )e^ X v ) + 0(|q/2 + k | 2 ) 



I 1 



O(|q/2 + k | 2 ) 



(53) 



where both and C v [Eq. (|43[) ] arc evaluated at 
k f« -q/2 and for fj, = {1, 2, 3, 4}, /x' = {2, 1,4, 3}. The 
Yafet part contribution I Y is not negligible for (/-process. 
The detailed derivation of I E and I Y (Appendix[F| bears 
some similarity with that of intravalley spin flip, and in- 
vokes the selection rules governing the opposite points of 
A star. 

We obtain the leading ^-process matrix element, 



M SiLA (k 1 ,ft;k 2 ,4) 



Dg S (K x 



iK y ), (54) 
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where 



D 



2P/E g:X [D so + 4A x k' D zz /A c } + O(k' 3 ) 



0.1 eV. 



(55) 



D zz = S,i + S u denotes the sum of dilation and uniax- 
ial deformation potential parameters [see Eq. (1E14[) and 
discussion after Eq. (|E16j) ] and D so ss 6.7 meV is a spin- 
dependent scattering constant [see Eq. (IF12I) ]. 

Squared matrix elements of the p-process are summa- 
rized in Table |Vj and Fig. [5] compares these analytical 
expressions with numerical results along typical K di- 
rections. The scattering-angle dependence is predicted 
correctly where the linear relation holds best when both 
states are at the vicinity of the valleys centers. From 
numerical results we note that O(k' ) 3 terms from Elliott 
interaction yield a negative 10% correction. 

Before concluding this part we briefly discuss spin-flips 
matrix elements of higher order. We only mention the 
general nature of these matrix elements without deriv- 
ing explicit forms. The reason is that the g-process has 
the weakest contribution to spin relaxation in unstrained 
bulk silicon. Among the higher-order matrix elements 
[last line of Eq. (|49|)]„ the 0(K 2 ) terms vanish by time 
reversal symmetry. This property can be proven by using 
Eq. p6[) and the matrix elements expansion to quadratic- 
K terms. On the other hand, the non-vanishing second- 
order term has a wavevector dependence of the type 
C(A")C(|q/2 + ko|). Such a symmetry allowed term is 
more important than 0(K 3 ) terms and it is governed by 
matrix elements of TA phonon modes. Its contribution to 
spin relaxation time from numerical calculations will be 
discussed in Sec. lVIII Here it is noted that compared with 
intravalley processes the relative contribution to spin re- 
laxation from high-order terms is larger. In addition, the 
signature of further bands (outside the X\ and X4 basis 
states) is larger in ^-process spin flips. 



D. Spin orientation dependence 

We relax the restriction of a fixed spin orientation 
(n||z) and explore this degree of freedom. The spin orien- 
tation dependence originates from the anisotropy of the 
conduction band. This anisotropy suggests that only the 
relative direction of spin orientation to a valley axis is 
important. Accordingly, we choose the +z valley and 
express n in terms of polar and azimuthal angles (6, (j>) 
with respect to the +z and +x directions. This choice 
is convenient since selection rules between A-point ba- 
sis states of the +z valley have been already derived in 
Sees. IVAWBl 

Derivation of spin-flip matrix elements when n | z 
relies on the specific form of |k, ft (ft) n )- These states 
satisfy the conditions of Eq. ([6]) and their forms can be 
found by solving the Hamiltonian matrix when written 
in terms of the new basis {|X) <g> | | (|) n )}- However, 
we can avoid this labor and find these states using a 



two-step procedure. First, the previously derived forms 
of |k, ft (ft) z ) [Eq. are re-expressed in terms of the 

new basis {|X)<g>| f (l)n)}- This change of basis amounts 
to a rotation of the spin coordinate by —9 about the axis 
u) = n x z. Specifically, the operator 



U = cxp 



(56) 



is applied on each of the four spinors [A Xi (k), B Xi (k)] T 
that Eq. (|27|) is comprised of.*^ 2 - In the second step, we 
find |k, ft (ft) n ) by forming linear combinations of the re- 
expressed states (still oriented along z) such that Eq. ^ 
is satisfied. Technical details of finding this superposi- 
tion are summarized in Appendix [D] The states can be 
expressed in the original {|X) ® | f (l)z)} basis using W . 

Having the forms of |k,ft (ft) n )j we repeat the proce- 
dure of Sec. IV Bl and derive the dominant intravalley spin- 
flip matrix elements. Scattering with acoustic phonon 
modes reads 



Mr ra (ki,ft n ; k 2 ,ft n ) 



iv/&cD' xy (fe£+ A (q) + g»£+ A (<J)) 



cos -{q x - iq y ) - 



sin — e 
2 



2 k; 



(q x + iq y ) 



. (57) 



where A={TAi, TA 2 , LA}. Scattering with optical 
phonon modes reads 

Mf^ftn; k 2 ,ft„) 
= - V /A c D op £" x (q) 



9 9 
cos 2 -{q x - iq y ) - sin 2 -e 211 ^^ + iq y ) 



(58) 



where x ={TOi, T0 2 , LO}. Table E] lists the squared 
matrix element expressions of all phonon modes using 
elastic continuum approximation for diamond crystal 
structures. These spin-flip expressions are specified for 
spin orientations along all inequivalent high-symmetry 
crystal directions. This chosen set of directions is im- 
portant for two reasons. First, the invoked clastic con- 
tinuum approximation is accurate along these directions. 
Second, the oriented spins in typical spin injection exper- 
iments point along these directions. Using the results in 
Table [V] we have shown in Fig. [5] that the analytical ex- 
pression for n||z (#=0) agrees with independent numeri- 
cal calculations. This agreement is also true (not shown) 
for n ft z cases in Table |Vj As will be explained in the 
next section, the numerical models automatically incor- 
porate time reversal and space-group symmetries of the 
crystal. The agreement between these independent cal- 
culation approaches manifest the robustness of the major 
terms we have kept in intravalley spin flips. 

Deriving the g-process spin-flip matrix elements is sim- 
ilar. We study scattering from z to — z valley for an arbi- 
trary direction of n. Repeating the procedure of Sec. IV Cl 
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the dominant spin-flip matrix element reads 
M gtLA {k u t n k 2 ,4 n ) = 



D 



cos 2 -(K x - iK y ) - sin 2 -e 2 ^(K x + iK y ) 



, (59) 



and it originates from scattering with LA phonon modes. 
Using this expression, Table [V] lists squared matrix ele- 
ments along high-symmetry crystal directions. As can 
be seen from the square brackets terms of Eqs. ([57|) - (|59l) . 
the ^-process shares the same angular dependence as in 
the intravalley case but with replacing q = ki — k 2 with 
K = ^(ki + k 2 ). The similar angular dependence is not 
surprising for two related reasons. First, the spin orien- 
tation only affects the electron states while the phonon 
properties play no role in setting the angular depcndance. 
Second, the electron states that we use in deriving in- 
travalley or ^-process matrix elements are all expanded 
around the same X point. The replacement of q with K 
can be understood by time reversal symmetry. 

Analysis of the spin relaxation time due to scattering 
within the +z valley and between ±z valleys will be pro- 
vided in Sec. I VII Cl for various directions of n. 



VI. INTERPLAY BETWEEN ANALYTICAL 
DERIVATION AND EPM 

In this work we compare our results with an empirical 
method in which the electronic states and phonon po- 
larization vectors are calculated, respectively, via empir- 
ical pscudopotcntial and adiabatic bond charge models 
(EPM and ABCM), 2 ^ The states and polarization vec- 
tors are used in calculating the electron-phonon interac- 
tion following a rigid-ion approximation^ The EPM and 
ABCM provide very accurate symmetry-related results 
and trends of contributions from high-order wavevector 
components. Given that a sufficiently large plane-wave 
basis is employed, then in addition to time reversal sym- 
metry these models capture the symmetries of the Bra- 
vais lattice and of the primitive cell. This ability is in- 
dependent of the specific chosen values of empirical pa- 
rameters (e.g, form factors of the pseudopotential). On 
the other hand, intensive numerical calculations do not 
automatically guarantee an accurate spin-flip matrix el- 
ement result. Our theory provides a clear insight to the 
identity of critical empirical parameters that are relevant 
for spin relaxation. In this section we will elaborate on 
fundamental aspects in understanding the application of 
EPM in spin-flip processes. 

The usual way of finding the energy band structure by 
adjusting the pseudopotential form factors is not suffi- 
cient for scattering problems. Specific derivative values 
of the pseudopotential at the first few reciprocal lattice 
vectors \d m V(k)/dk m at k = g n ] are additional neces- 
sary conditions. In momentum scattering, a correctly 
interpolated pseudopotential is capable of reproducing 
the energy shifts of the conduction band in response to 



applied stress. Spin scattering is more than momentum 
scattering in the sense that the leading-order matrix ele- 
ment is of higher order in the wavevector (e.g., intraval- 
ley and g-process scattering in silicon). One consequence 
is that other deformation potential constants may come 
into play (e.g., D' xy as was shown in the previous sec- 
tion). In fact, our spin-dependent EPM is matched not 
only with energy band structure but also with different 
deform ation potential quantities (£ i + a, b, d, £2 and £| 
in Ref. (lO^ L 

The analytical derivation shows that the two-band de- 
generacy at the X point plays an important role in silicon 
due to its proximity to the valley center. This proximity 
results in Ac /E g X -C 1 which allowed us to discard inter- 
action terms of the type r i rj9V+ / drk [Eq. (|42f|) ] . Consid- 
ering the Fourier transform of the pseudopotential, this 
simplification means that intravalley spin flips are only 
sensitive to values of V(g n ) and dV(g n )/dk while high- 
order derivatives can be discarded. However, silicon is a 
specific case. In germanium, for example, intravalley spin 
flips have a cubic dependence on the wavevector ji 2 ^ and 
they may require additional information on d 2 V (g n ) / dk 2 
and d 3 V(g n )/dk 3 . In the framework of deformation po- 
tential theory, this amounts to expanding the strained 
crystal potential with quadratic or higher-order strain- 
tensor components. Notably, a comprehensive experi- 
mental analysis of the intravalley spin relaxation time in 
germanium may provide new information on its crystal 
potential. 

The /-process spin-flip matrix elements in silicon de- 
pend on various parts of the pseudopotential curve. An 
a priori and independent determination of V(|k|) is dif- 
ficult. To empirically interpolate the pseudopotential 
curve from momentum and spin relaxation experiments, 
one has to know the relaxation times dependencies on 
temperature, electric field, stress and related 'knobs'. 
This information can resolve the values of individual scat- 
tering constants rather than their combined effect. In ad- 
dition, knowledge of the intravalley scattering parameters 
allows one to reduce the uncertainties in interpolating the 
pscudopotcntial. 

All in all, the kp method and group theory provide an 
unambiguous guidance in relating EPM parameters with 
different experimental measurements and in determining 
the relevant parts for spin relaxation. These consider- 
ations rationalize the investigation of spin-flip problems 
by joining analytical and numerical approaches. 



VII. SPIN LIFETIME 

The spin relaxation time is an experimentally acces- 
sible quantity. With a specific electron distribution J 7 , 
detailed expressions for spin-flip matrix elements can in 
principle provide a transparent physical picture of the 
spin relaxation under a variety of conditions. The spin 
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relaxation rate has the form. 




where (k 2 |Hpp|ki) denotes the expression of Eq. (UJ) with 
opposite spins Si and S2. The material volume Na 3 /A is 
chosen as the unit volume. +(— ) corresponds to phoiion 
emission (absorption). The average over ki represents 
BF/dE^ weighted integration over ki, which is exact 
at the limit of infinitesimal spin-dependent chemical po- 
tential splitting. The prcfactor of Air/K instead of 2ir/K 
denotes the fact that the net number of spin-polarized 
electrons (N^-N^) changes by two with each spin flip. 

A few applications will be shown in this section. Af- 
ter briefly describing our numerical integration effort, we 
present the commonly used r s f under a normal condi- 
tion both analytically and numerically. Here 'normal' 
stands for a non-degenerate bulk silicon without strong 
fields (i.e., Boltzmann distribution Jmb for electrons), 
and with spin orientation along the valley axis. The av- 
erage over ki in Eq. (pO)) then becomes J-"mb weighted 
integration over ki. This weighted integration is valid 
thanks to the relation cU~mb / dE oc Jmb- In the last 
part of this section, some essential relations between the 
spin orientation and r s f are derived. 

A. Numerical integrations and approximations 

We have performed numerical integrations of Eq. ([50)1 
at different levels of approximation. These calculations 
are presented in decreasing order of their computation 
time. 'EPM+ABCM' denotes the full numerical results 
from EPM and ABCM program codes^ The calculated 
electron states and phonon polarization vectors are then 
incorporated into a rigid-ion model following the proce- 
dure in Ref. [H[. 'kp+ABCM' replaces the EPM results 
with kp energy band [Eq. (|C4[) ]. It employs the general 
analytical form of |M sf | from Eqs. gS]), (@TJ) and ([51)) . 
'cllip+u;(g)' further employs a spheroidal energy disper- 
sion and replaces the phonon frequency from ABCM with 
lo = gi'TA/LA f° r intravalley and g-process scattering with 
acoustic phonon modes (wta/la are the phonon veloci- 
ties) . For intravalley scattering with long- wavelength op- 
tical phonon modes, it replaces the phonon energy with a 
constant fkjj op —63.5 mcV. Finally, it employs the elastic 
continuum approximation for phonon polarizations (Ta- 
ble H with n|| [001]). 

Numerical integrations of Eq. (|60[) were performed us- 
ing a grid spacing of 0.01 x 27r/a in k space. This 
grid leads to converged results for intermediate and high 
temperatures (T > 50 A'). We take advantage of the 
eight-fold symmetry of the A axis valleys and reduce 
the intensity of computation whenever possible. The 
edges of the irreducible wedge are weighted to pre- 
vent overlap with their neighbors. Also helpful is the 
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FIG. 7: (Color online) Intravalley r s f (T) induced by scattering 
with TA, LA and optical (OP) phonon modes. The relaxation 
time is contributed from all valleys and the spin orientation 
is set along one of the valley axes (A axis). Analytical curves 
of the acoustic modes follow Eq. (|61|) . Other curves refer to 
numerical approaches with calculation intensities that depend 
on details of the band structure, phonon energy, and polar- 
ization vectors (see text). Due to the crossing of TO and LO 
dispersion curves, we do not separate their contributions in 
the ABCM calculation. 



strict equality j r MB(ki)?i(q) = J r MB(k2)[^(q) + 1] when 
E(k.2) — E(ki) = hu>(q). Typical execution times are as 
follows: tens of seconds for 'ellip+cj(q , )' ) tens of minutes 
for 'kp+ABCM' and a few days (with 64 CPU cores) for 
'EPM+ABCM'. The bottleneck of the 'EPM+ABCM' 
computation speed lies in calculation of individual matrix 
elements. This calculation involves the product of states 
and interactions written, respectively, in vector and ma- 
trix forms with a basis of hundreds of plane waves 



B. Fixed spin orientation along the valley axis 



We separately study the spin relaxation time due to 
intravalley, g-process, and /-process spin flips. Figure [7] 
shows results of the intravalley spin relaxation time by 
integration of Eq. (|60| using the aforementioned nu- 
merical procedures. The figure also includes analytical 
curves for scattering with acoustic phonon modes (left 
and middle panels). These analytical integrations are 
carried out by employing an elastic scattering approxi- 
mation and a high temperature limit for phonon popu- 
lation n(q) sa k^T /huj(q). Together with the mentioned 
simplifications in numerical integrations, these common 
practices allow us to accurately calculate the relaxation 
time by employing the relevant spin-flip matrix elements 
in Table|V](|A4 n ^ A | 2 with n|| [001] and n||[100]). The 
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integrated spin relaxation rate reads 
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where and r s L f A j are, respectively, the intravalley spin 
relaxation times due to scattering with long-wavelength 
TA and LA phonon modes. Table I VIII lists the values 
of all parameters in the above expression. To enable an 
accurate analytical integration in comparison with the 
full numerical integration, the explicit dependence on the 
band structure anisotropy has been considered*^ This 
anisotropy is expressed in Eq. (|61|) via r — m t /mi which 
denotes the ratio between the longitudinal and transverse 
effective masses of the electron (with respect to the valley 
axis). As seen from the middle and left panels of Fig. [7j 
the analytical integrations match very well with the most 
detailed numerical integration. These figures also show 
that l/r sf (T) of the 'EPM+ABCM' numerical results de- 
crease slightly faster than T~ 5 / 2 for both scattering with 
TA and LA phonon modes. It indicates the dependence 
of |Af s f | on higher-order wavevector components when |q| 
gradually increases j 21 - The figure shows that only minor 
changes are introduced in all intravalley processes when 
replacing the detailed conduction band structure with 
spheroid dispersion and the numerical phonon data with 
the analytical approximation. This behavior supports 
the validity of the invoked approximations, including the 
use of an elastic continuum approximation. The stronger 
deviation of the 'EPM+ABCM' curve in the case of opti- 
cal phonon modes shows the effect of higher-order matrix 
element terms. Although the contribution to spin relax- 
ation from LA and optical phonon modes is negligible 
at low temperatures (compared with TA modes), their 
effect should be considered at room temperature (espe- 
cially the optical modes). 

We study the spin relaxation time due to g-process 
spin flips in a similar way. Results of the analyti- 
cal and numerical integrations are presented in Fig. [3] 
The analytical relaxation rate due to scattering with LA 
phonon modes (dominant effect) is reached by integrat- 
ing Eq. (|60[) with spin-flip matrix elements taken from 
Table[V](for the case of n||[001] and n||[100]), 



V2D 2 gs m t mjEI 



K 2 (E$/2k B T) 



'sf,9 



16ttS?1 4 p VkBTsinh(El/2k B T) 



(62) 



Values of the scattering constant, phonon energy and ef- 
fective masses are listed in Table [VTTl Ki is the i th order 



TABLE VII: Parameter values in Eqs. J6TJ and (J62J) 
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FIG. 8: (Color online) g- process r B {(T) from LA phonon at 
different levels of approximation (see text). The 'TA' curve 
is calculated only by 'EPM+ABCM'. The relaxation time is 
contributed from all valleys and the spin orientation is set 
along one of the valley axes. In 'ellip+t<j(g)', the phonon ve- 
locity is the same as that of a long-wavelength LA phonon 
mode. In the analytical curve the phonon energy is constant 
(21 meV). 



modified Bessel function of the second kind. Figure|8]also 
shows a full numerical curve due to scattering with TA 
phonon modes. The relative weight of the TA part clearly 
increases with temperature. Also can be seen at the high 
temperature end, is a relatively large deviation between 
the analytical and 'ellip+w(g)' curves. This deviation is 
caused by using a constant phonon energy instead of a 
linear dispersion relation. The energy difference affects 
the phonon population. 

Finally, we study the spin relaxation due to /-process 
spin flips. Analytical and 'EPM+ABCM' numerical in- 
tegrations of Eq. ([60)) are presented in Fig. [9] for each 
of the phonon modes. In the analytical integration, we 
have used the approximations of a spheroidal energy dis- 
persion in the conduction band, and of wavevector inde- 
pendent spin-flip matrix elements and phonon energies 
(Sec. HVj) . The analytical integration becomes relatively 
simple for each of the nonvanishing modes and for spin 
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orientation along the valley axis, 

4 = ^Ud1 K ^^f ] • (63) 

Values of the scattering constants and phonon energies 
(D-£ i and E?) are listed in Table Hill where stands 
for The value of Ai is obtained by summing |M s f | 2 

over 12 pairs of valleys. Following the results of Sec. II V Al 
we get that A± = 8 and A 2 /z = 16. In calculating the 
values of Ai we have used the facts that for £i, there is 
no coupling between valleys of the x-y plane whereas the 
coupling is D|, between any of the remaining eight pairs 
of valleys (assuming that n|| z). For S 2 /3, the coupling 
is |(1 + i)Dz 2/3 \ 2 between each of the 4 pairs of the x-y 
plane and £)|, between each of the remaining 8 pairs. 

Figure [9] shows that the dominant scattering is with 
phonon modes of the LA and TAy (Ei and £3 symme- 
tries; see Table [TTT]) . Their temperature trends are cor- 
rectly predicted by the wavevector-independent analyti- 
cal analysis. The same applies for scattering with phonon 
modes of the TOj_ and TOy. After a qualitative anal- 
ysis we find that the big difference between analytical 
and numerical results in the LO phonon case is caused 
by quadratic- wavevector terms. This wavevector depen- 
dence is further complicated by even higher-order terms 
when the electron states are further away from valley cen- 
ters. For the TAj_ phonon case, we found in Sec. IIVI that 
there is no wavevector independent term. The numerical 
result is attributed to linear terms and it is non-negligible 
due to the higher phonon population of this mode [low- 
est energy; see Fig. |3ja)] . In general, wavevector depen- 
dent terms contribute at all modes when the temperature 
increases. The dependence on the electron wavevectors 
(ki^) can be similarly analyzed between decomposed kp 
basis states (as in the intravalley and <?-process cases). 
We do not make an explicit derivation of these terms 
since the focus is on the leading-order term contribution 
(wavevector independent in the /-process case). 

C. Spin orientation anisotropy coupled with 
symmetry breaking mechanisms 

The spin relaxation time is in general a function of spin 
orientation [r s f(n)]. It can be obtained by integration 
with the general matrix element expressions derived in 
Eqs. (Tr?a) - (fT7f]) and We can firstly calculate 

r s f(n) from one valley (intravalley) or a pair of valleys 
(intervallcy) . The total r s f(n) is a summation from all 
of the involved valleys, whose individual r s f(n) can be 
related to the calculated one by a proper coordination 
rotation. 

We first discuss the integrated effect of changing n on 
intravalley and ^-process spin flips. Figure [TOl shows the 
temperature dependence of r s f (n) due to intravalley scat- 
tering of electrons in the z-valley with TA, LA, and op- 
tical modes. It also shows r s f(n) due to the dominant 
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FIG. 9: (Color online) Analytical (solid lines) and numerical 
(-x- lines) results of r B f(T) in /-process spin flips where each 
phonon mode is represented by a different color. The relax- 
ation time is from all valleys and the spin orientation is set 
along one of the valley axes. _L and || are taken with respect 
to the cross product of valley centers (direction of ko,i x ko,2). 
The analytical contribution from TAx modes vanishes (at the 
zeroth order). 



intervalley scattering between the ±z valleys (g- process). 
These results were calculated using the 'kp+ABCM' inte- 
gration procedure with matrix elements taken from Ta- 
ble [V] Each panel shows results of n along all of the 
inequivalent high-symmetry crystallographic directions. 
The spin lifetime in a given valley increases with decreas- 
ing the projection of the spin orientation on the valley 
axis. This effect implies that suppression of the spin re- 
laxation in one valley is compensated by enhanced re- 
laxation at perpendicular valleys. As a result, the spin 
lifetime due to intravalley and g-process scattering from 
all valleys is expected to have a diminished dependence 
on the spin orientation. Nonetheless, Fig. [TU]can be seen 
as a simplified example of how different spin orientations, 
coupled with a symmetry-breaking mechanism, can lead 
to different experimentally measurable quantities. For 
example, stress or electrical fields can selectively change 
the electron distribution of different valleys leading to 
a significant electron population only at valleys along a 
certain axis. In that case, the total spin lifetime will 
present a pronounced degree of anisotropy when chang- 
ing the spin orientation. 

For /-process spin flips, changing the spin orientation 
results in slightly more involved relations between spin re- 
laxation times of different pairs of valleys. The wavevec- 
tor integration of Eq. ([ST)]) is not affected by the spin ori- 
entation (M^ f are wavevector independent). Therefore, 
one only needs to obtain the values of A, in Eq. (j63| for 
each of the £j symmetries. In this application, the 12 
pairs of involved valleys can be divided into 6 groups. 
Each group consists of two pairs related by space inver- 
sion operation and they always have the same value of 
|M S4 (n)| = (k 2 ,4 n |.ff Ei |ki,ft n ). For ±x O ±y pairs, we 
use Eqs. (|17bj) . (|17d|) and (fT7fj) . and denote them with 
Mo^^n). Using crystal symmetry, the matrix elements 
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FIG. 10: (Color online) r a t(T, n) of electrons in the z valley 
(intravalley) and ±z valleys (g-process) for spin orientation n 
along high symmetry directions. The curves of n|| [001] and 
[110] overlap each other since integrations of their respective 
spin-flip matrix elements over the azimuthal angle of q yield 
the same result. In this sense, only the projection h z on the 
valley axis is relevant. 



of all other pairs relate to Mq^, (n) by, 



±x <-» ±z 
±x 
±y «-» ±z 



M Sj (n 



M ,s j (n x , -n y ,n z ), 
M ^ i (n x ,n z ,n y ), 
M a ^ t (n x ,-n z ,-n y ), 
M ^(n z ,n y ,n x ), 
M ^ t (-n z ,n y ,-n x ). (64) 



After summing all of the contributions, one finds that 
A\ = 8 and A 2 /s = 16 [the same as those obtained 
in Eq. (|63[) ] with the direction of n first expressed in 
terms of (#, (f>) and then substituted into expressions of 
Mq,^ ■ Therefore, when all valleys are equally populated 
then the total r s f is invariant of n (to the leading order). 
On the other hand, when symmetry breaking effects are 
introduced then valley repopulation brings in a depen- 
dence of the spin relaxation time on the spin orientation. 
For example, consider an internally strained structure in 
which the /-process scattering with ±2 valleys is sup- 
pressed (i.e., the ±x and ±y valleys are equivalent and 
have sufficiently lower energy). In such a structure, the 
above analysis reveals that if the spin is oriented along 
z then one should assign in Eq. (|63p values of A\ = 
and ^2,3 = 8. Similarly, for orientation along x or y one 
should assign A\ = 4 and A 2 .z = 4. Using the param- 
eters of this stressed configuration, the /-process spin 
relaxation time is estimated to be shortened by ^50% 
when changing the spin orientation from the z-axis to 
the perpendicular plane. 



VIII. SUMMARY AND OUTLOOK 

We have presented a comprehensive analysis of all 
phonon-induced spin relaxation processes in bulk sili- 
con. The applied temperature and doping regime of this 
mechanism has been identified, among other mechanisms 
(Sec. |IT]). I n decreasing order of contributions to spin re- 
laxation, detailed expressions of /-process, intravalley 
and (/-process matrix elements have been derived and 
their dependence on the spin orientation are unveiled. 
We have elaborated on the wavevector dependence and 
symmetry properties of each spin-flip process. In analogy 
to Herring and Vogt theory on momentum relaxation in 
siliconfiSi this work unravels the magnitudes and sym- 
metries of all phonon-induced spin relaxation processes 
in silicon. 

In studying the /-process spin flips, double group selec- 
tion rules are used to obtain the wavevector-independent 
(leading order) matrix elements. Spin orientation depen- 
dent spin-flip (and spin-conserving) matrix elements were 
expressed in terms of scattering constants D^ iS (and also 
-Dsim), for electron-phonon interaction with symme- 
try. 

Intravalley spin flips were studied by using a combina- 
tion of single group theory, kp perturbation method and 
rigid-ion model. The spin-dependent coupling between 
the expanded basis functions and symmetrized interac- 
tions is formulated via selection rules. This approach 
allowed us to derive the leading order intravalley matrix 
elements [Eqs. ([57]) and ([55]) ] and to resolve their exact 
dependence on the phonon polarization (£), its wavevec- 
tor (q), and on the spin orientation with respect to the 
valley axis [n(0, <fi)]. By incorporating the diamond crys- 
tal structure into an elastic continuum model (expressing 
£ in terms of q) , we have derived appealing forms of these 
matrix elements (TablelVf. Finally, the analysis identifies 
the important band structure parameters. The r\ parame- 
ter is a measure of the wavevector independent spin-orbit 
coupling between conduction and valence states at the X 
point (Table IVI[) . The deformation potential parameter 
due to scattering with long-wavelength acoustic (optical) 
phonon modes is D op (D' xy ), and it corresponds to in- 
terband coupling between the lowest pair of conduction 
bands [Eq. (gBJ) and (|4"8"])]. This coupling is brought by 
the proximity of the valley to the two-band degeneracy 
at the X point. The Ac parameter denotes the energy 
gap between the lowest pair of conduction bands at the 
valley center. 

A complete picture of the relation between intravalley 
spin-conserving and spin-flip processes has been provided 
(together with the detailed derivations in Appendix [E] 
and [F]) . This comparison reveals important physical as- 
pects that are being overlooked when relating spin and 
momentum relaxation times via the shift of the g-factor 
(conventional approach in quantifying the spin lifetime 
due to the Elliott-Yafet relaxation mechanism). 

g-process spin flips are studied in a similar way to the 
intravalley case. In spite of the opposite valley positions 
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of an electron before and after scattering, Kramers conju- 
gation relation allows us to expand the electronic states 
by basis functions of the same X point. One result of 
this relation is that the matrix elements depend on the 
average between the initial and final electron wavevec- 
tors [K = ^(ki + IC2)] rather than on their difference 
(q = ki — k2). The spin orientation dependent ma- 
trix element of the g-process is provided in Eq. (|59|) . 
The involved scattering constant has a large contribu- 
tion from dilation and uniaxial deformation potential 
constants [Eq. (|55[) ]. Comparing the derived matrix el- 
ements with respective results of independent numerical 
calculations shows that our analytical approach provides 
accurate spin-flip amplitudes at all scattering angles for 
both intravalley and ^-process cases (Figs. [5] and [5]). 

Our analysis provides insights into which parts of the 
interaction (spin-independent 'Elliott' or spin orbit cou- 
pling 'Yafct') dominate the phonon- induced spin relax- 
ation in silicon. In silicon, the sum of Elliott and Yafct 
contributions vanishes at the zero and first order of in- 
travalley scattering with acoustic phonon modest At 
the leading order of this scattering (quadratic-in-q), we 
have shown that the Elliott part dominates the spin re- 
laxation. The Elliott contribution is also shown to domi- 
nate the spin relaxation due to intravalley scattering with 
optical phonon modes and g-process intcrvalley scatter- 
ing. The latter two process are, respectively, linear in 
q and K. Yafet contributions, however, cannot be com- 
pletely ignored in g-process scattering where the defor- 
mation potential constant is affected by the spin-orbit 
coupling [Eq. (|F12[) ]. In silicon, Elliott and Yafet contri- 
butions are comparable only in the /-process. 

We have derived the spin lifetime due to each of the 
spin-flip processes by integrating its leading matrix ele- 
ments. Analytical forms are given for intravalley scat- 
tering with acoustic phonon modes [Eq. (|61[) ] and for 
both types of intervalley scattering [Eqs. (|62|) and (|63|)]. 
Comparison of these results with numerical calculations 
at different levels of approximation show good agree- 
ments (Figs. [3 [8] and [9]). The analysis also identifies 
the phonon modes which lead to strongest spin relax- 
ation. The /-process is led by scattering with £1 and 
S3 symmetry (LA and TAy) phonon. Intravalley spin 
relaxation is led by scattering with TA phonon. Intraval- 
ley contributions from scattering with LA and optical 
phonon become comparable to the TA's at room temper- 
atures, g-process is led by scattering with LA phonon. 
We have also considered the secondary contribution to g- 
process spin relaxation from scattering with TA phonon 
(quadratic-wavevector dependence but a larger phonon 
population). 



Outlook 

Results of this work shed light on new research di- 
rections in group IV spintronics. By having a thorough 
understanding of the underlying physics, one can devise 



a means to enhance the spin lifetime of room tempera- 
ture silicon spintronic devices. Quenching of the domi- 
nant /-process by certain stress configurations is the first 
step in this direction! 15 i 22 In this case, only the valleys 
along one crystallographic axis are practically populated 
with electrons. To further suppress the remaining re- 
laxation processes (intravalley and g- process spin- flips), 
one can make use of the slower relaxation when the spin 
orientation is perpendicular to this crystallographic axis 

(Fig. on]). 

In addition to quenching the /-process and further 
optimization by spin orientation, one can also impose 
geometrical constraints on the transport. For example, 
promising candidates seem to be stressed silicon wires 
with a cross-section area that is large enough to prevent 
detrimental surface effects but is small enough to restrict 
the phase-space for scattering. Having the wire axis par- 
allel to the axis of populated valleys would allow one to 
achieve significantly longer spin lifetimes. This fact can 
be seen from the detailed intravalley and g-process ma- 
trix elements: while in these structures forward and back- 
ward scattering with respect to the axis of the wire (and 
populated valleys) dominate the transport, these types of 
scattering would not be accompanied by spin flips (e.g., 
for z valley electrons, assign q x m q y as in Table fV|) . 
This example shows the insights one can gain from un- 
derstanding the symmetries of the matrix elements rather 
than only having a knowledge of the integrated effect. 

The complete set of matrix elements is also instrumen- 
tal in calculating the spin relaxation in the presence of 
large electric fields. Due to the mass anisotropy, the val- 
ley population depends on the direction of the field. In 
addition, the field can lead to a large departure of the 
electron distribution from equilibrium conditionsJ^ As 
a result of these effects, intervalley processes are enabled 
already at low temperatures and certain scattering pro- 
cesses are enhanced^ Using the dependence of the relax- 
ation on spin orientation and scattering directions, one 
can accurately model and understand the spin relaxation 
in these conditions. 

The presented theory identifies a handful of scatter- 
ing and band structure constants which have not been 
experimentally determined yet. Evidently, the most im- 
portant constants are Ds 1(3)S of the /-process scattering 
with phonon modes along the LA (TAy) branches, and 
Ax which denotes the spin-orbit coupling between con- 
duction and valence states at the X point. In the absence 
of experimental data, we have used the empirical pseu- 
dopotcntial method to calculate their values (Tables IIIII 
and IVI[) . To determine these constants experimentally 
one should resolve various contributions to the measured 
spin lifetime. Intravalley, g and /-processes have different 
dependencies on the wavevector components, phonon po- 
larization and spin orientation. In addition, the energies 
of the respective phonon modes are different. As a result, 
the measured spin lifetime of each of these processes has a 
unique dependence on temperature [Eqs. (|6I D - (l63p ]. and 
it can be clearly resolved by application of a symmetry 
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breaking mechanism (e.g., stress or electric fields as men- 
tioned before). 

Finally, the theoretical approaches presented in this 
paper can be used to study the spin relaxation of ma- 
terials with different symmetry groups and consequently 
different wavevector-order analysis (e.g., germanium and 
graphene with respective utilization of the space groups 
at the L and K points of their Brillouin zones). Results 
of such a study provide a clear picture of preferred scat- 
tering angles, spin orientation and dominant spin relax- 
ation mechanisms. As in the case of silicon, having this 
information provides guidance in tailoring the spin relax- 
ation by application of stress, external fields or geomet- 
rical constraints. When such external influences become 
too large, one can repeat the steps of the presented pro- 
cedure after adding the external perturbation explicitly 
in the Hamiltonian and interaction terms. 

This work is supported by DOD/AF/AFOSR FA9550- 
09-1-0493 and by NSF ECCS-0824075. 

Appendix A: detailed application of double group 
theory for /-process matrix elements 

In this appendix, we intend to express interaction ma- 
trix elements (k2,S2|i?s i |ki,Si) in terms of Af^ inde- 
pendent constants. First, by the general time reversal 
and space inversion symmetries we can connect different 
matrix elements of each phoiion mode [Eq. (|13[)]. 

For Si mode, the {pz\i~) operation equates the spin- 
flip matrix element to negative of itself (seen from the 
character of Si and the IR matrix of _Da 6 = -Dai x -D1/2 
in Table IIII[) . We show this example explicitly, 







k,1r») - 


-> e 4 x (-i)\k,f z ) 




■» e 1 x z|k,4J. 2 ), 



remembering that basis states in Z?a 6 are |k,l"|-z (JJ- Z )). 
(p z \ T ) an d other symmetry operations do not provide 
constraints on the spin-conserving matrix element. Since 
A/si = 2 [Eq. ([IT])], there are two real constants D^ 1 , m 
and Ds 1>s . The physical significance of these two con- 
stants will become clear later. They are defined such 
that 

(k 2> ft, l-ffsjki.ft,) = -D Sl ,m + iD EliS . (Al) 

For £2, spin-conserving matrix element vanishes (can 
be seen, for example, by applying the Ta z operation). 
The spin-flip matrix element is 

(k 2 ,^ |# S2 |ki,lU T = -(-ki.-Mtfal-ka.ft,) 

ip ^ T} i(ka,-fr, |^ 2 |ki,^), 

by time reversal and (p xy \t) sequentially. Combined with 
Eq. (|I3b|) . we are led to 

(k 2 ,4z|ff£ 2 |ki,lU = £>s 2 -;£>£ 2 , (A2) 



where Ds 2 is the independent real constant. 

Applying the same operations for S3 and noting 
that the only difference from £2 mode is the sign of 
Xt. 3 {Px V \t), we get 

(k 2 ,4 2 |#£ 3 |ki,^) = D^+iD^, (A3) 

where Z?e 3 is the independent real constant. As shown 
in Eq. (jlip . the £4 phonon symmetry does not couple 
electrons of any spin species (A<£ 4 = 0). 



Appendix B: Group G\ 2 

Symmetry operations of the G\ 2 group are listed in 
Table IVnTJ The table refers to the X point in the z di- 
rection. Tables of X points in the x and y directions 
are derived by cyclic permutations. From Table [Villi one 
can study how vectors and axial- vectors are transformed 
under the 32 group elements. Character table of the G§ 2 
group is listed in Table IIXI From the character table one 
can decompose the direct products of {X\, X4} into di- 
rect sums of IRs [Eq. ([T5])]. One can also construct all 
sorts of direct product rules from IRs that present com- 
ponents of vectors, axial-vectors or (anti) symmetric po- 
tentials [Eq. ([20]) ]. 



TABLE VIII: List of operations of the G§ 2 group. Notations 
follow Ref[H. 



Class symbol 


Operations 


Ci 


(e|0) 


c 2 


(S 2x \0), (S 2 y\0), (S 2x \t x y), (S 2y \t xy ) 


C z 


S 2z 


C 4 


(5 2 xy\T),(S 2 xy\r + t x y) 


c 5 




c 6 


{PAT), (pz\T + txy) 


c 7 


(Pv\t), {px\r),{p y \r + t xy ), (px\r + txy) 


C 8 


(i\r), (i\r + t xy ) 


c 9 


(Pxy\0),{pxy\0) 


C10 


(ff4*|0), ((T4z\0), (<J4z\t xy ), {V4z\txy) 


Cn 


(pxy\txy) ,(pxy\txy) 


C12 


(S 2 xy\r),(S 2x y\r + t x y) 


Cia 


(5 2z \t X y) 


C14 


(e\txy) 



Appendix C: k ■ p Hamiltonian parameters and the 
Partitioning 

Comparing directly the matrices form of Eq. (|24|) and 
the Hamiltonian Hi in Eq. (f2"3"|) . the four independent 
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TABLE IX: Character table of the G| 2 group 



hi 


1 


4 


1 


2 


4 


2 


4 


2 


2 


4 


2 


2 


1 


1 




Ci C2 C3 C4 C5 C'a CV C8 Cg C10 Cn C12 


C13 C14 


Mi 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


M 2 


1 


1 


1 


-1 


-1 


1 


1 


1 


-1 


-1 


-1 


-1 


1 


1 


M 3 


1 


-1 


1 


-1 


1 


1 


-1 


1 


-1 


1 


-1 


-1 


1 


1 


M 4 


1 


-1 


1 


1 


-1 


1 


-1 


1 


1 


-1 


1 


1 


1 


1 


M 5 


2 





-2 








2 





-2 














-2 


2 


M[ 


1 


1 


1 


1 


1 


-1 


-1 


-1 


-1 


-1 


-1 


1 


1 


1 


M' 2 


1 


1 


1 


-1 


-1 


-1 


-1 


-1 


1 


1 


1 


-1 


1 


1 


M' z 


1 


-1 


1 


-1 


1 


-1 


1 


-1 


1 


-1 


1 


-1 


1 


1 


Mi 


1 


-1 


1 


1 


-1 


-1 


1 


-1 


-1 


1 


-1 


1 


1 


1 


M 5 


2 





-2 








-2 





2 














-2 


2 


Xi 


2 





2 

















2 





-2 





-2 


-2 


x 2 


2 





2 

















-2 





2 





-2 


-2 


x 3 


2 





-2 


2 























-2 


2 


-2 


x 4 


2 





-2 


-2 























2 


2 


-2 



integral constants readily follow, 

hk = (x 2 ' \p*\x?), 

m P = h(x£] Py \X?}, 
4m%c 2 A x = ih(x% |(W x p) y \X? 
4mlc 2 a = ~ih 2 (x% \V iV\ Xf' 



(Cla) 
(Clb) 
(Clc) 
(Cld) 



These constant can also be expressed as other equivalent 
integrals (e.g., hk = - (X{ \p z \ X{)). 

To analytically diagonalize the Hamiltonian matrix in 
Eq. (|23|) we note that the X-point energy gap, E g , x sa 
4.3 eV, is significantly larger than other energy scales. 
Therefore, we use degenerate second-order perturbation 
theory and lump the valence band effect onto the con- 
duction band (Lowdin partitioning) We get a reduced 
4x4 matrix, 



Hrr = H r ., 



E 9,X 



(C2) 



whose four spin-dependent basis states are 



x L = x 1 + 



-Xi 



(C3) 



Xi = [\x?,t),\x?,i),\XiA),\xhl)] T and x * = 

t>, \Xl I), \Xl t), \Xf, i)] T . Higher-order pertur- 
bation does not bring dominant terms up to quadratic k 
dependence. It has been explicitly checked. The eigen- 
values of H cc that pertain to the energies of the upper 



and lower conduction bands read 

h% 2 h 2 (k 2 > >» 2 



E ± (k) 

Zttiq 

AE c (k) = 2 x 



2m t 



k 2 )^ AE c (k) 



(C4) 
(C5) 



'h 2 k n k',\ 2 fh 2 k x k v \ 2 4A c A x ak' z 
1 'Wiki + k 2 )- 



m 



in,- 



E, 



9,x 



where k' z = k z — 2?t / 'a. The energy gap between the con- 
duction bands is AE c (k). Values of the parameters are 
provided in Table I VII The 77 related terrr>2£ as well as the 
ak' z related term in Eq. (|C4[) are the leading spin-orbit 
effect on the energy. Along the Z-symmetry axis connect- 
ing the X and W points \k' z = k y = or k' z = k x = 0; 
see Fig. [3Ja)] as well as the part of A axis very close to 
the X point, the splitting between the lower and upper 
conduction bands is induced by the spin-orbit coupling. 
This feature corresponds to the celebrated spin hot-spot 
at the edge of the Brillouin zonei 20 ' 21 Figures@|b) and (c) 
show the energy dispersion of the two conduction bands 
along the Z axis. 



Appendix D: Spin alignment 

In this appendix, we present the routine for linearly 
combining doubly degenerate eigenvectors such that the 
resulting eigenvectors spins are aligned along a desired 
direction of n [Eq. ©]. This routine is for Hamiltonians 
that include spin-orbit coupling and when the crystals 
have a space inversion symmetry. Particularly, it is tai- 
lored for basis states which go back to themselves after 
sequential space inversion and time reversal operations. 
The chosen X-point basis as well as general plain-wave 
basis belong to this category. Otherwise the routine can 
be readily modified and made applicable. 

In general, we express the double degenerate eigen- 
vectors in basis |X) <X> | j" (|) n ), where |X) is the spin- 
independent part (Sees. IV Al and IV D[) . <r • n in Eq. © 
is written as [1,0;0,-1] in this basis. Suppose one of 
the double degenerate eigenvectors is c a (k), and it has 
2m elements, where m is the number of spin indepen- 
dent basis states (m = 4 in our k • p Hamiltonian). 
c 2i-i(2i) (k) are coefficients of pure spin up (down) basis 
states. Then, from the general consideration of time re- 
versal and space inversion symmetries of the Hamiltonian 
and the basis states, we know that components of the 
other eigenvector c h (k) can be written as c b 2i = (cfj-i)* 
and C2i_i = — (c^i)* ■ This property satisfies the hrst 
equality of Eq. (|6]). To satisfy the other spin alignment 
definition in Eq. ((6]), we write a normalized linear com- 
bination 



,t _ 



such that, 



= (c° 



I m <g) 



wc 



')/vT 



(Dl) 



1 

-1 



2^ =0, 
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(4-1)*: 4-1 -(4)*. ^d/ m i 



fr 



whore c 

dimensional identity matrix. We are led to 

-d-wb + w 2 d* = 0, 

where 

TCI TCI 

&=E(i c 2i-ii 2 -Ki 2 )^=E( c 2i-i4)- 

i=l i=l 

Thus, we get the combination parameter in Eq. (|D1[) 



w = (6 - y / & 2 +4|d| 2 )/(2d*). 



(D2) 



We further have u> ss — d/6 if 6 3> d. In the case of EPM 
states, this general procedure requires to replace the |X) 
basis (m=4) with a plane-wave basis (where typically 
m>100). 



Appendix E: Intravalley momentum scattering and 
deformation potential 

In this appendix, we derive results of intravalley mo- 
mentum scattering by selection rules with k ■ p Hamilto- 
nian eigenstates. The procedure paves the way for ana- 
lyzing more involved intravalley (Sec. IV Bl) and (/-process 
(Sec. IV CP spin-flip scattering in a similar approach. 

The leading terms in momentum scattering matrix el- 
ements depend linearly on the phonon wavevector. We 
show it by a wavevector-order analysis. 

Zeroth-order: At = k 2 , the in-phase atomic vi- 
bration does not depend on lattice sites. Thus, the 
phonon-induced interaction in Eq. reduces to dis- 
placement of the entire crystal, V r V+(r — Rj) = 
VV cr ystai- Based on the relation 



VV, 



crystal 



[ P> H]/h, 



(El) 



the coupling of the in-phase part between spin-degenerate 
eigenstates of H vanishes, 

Mi(k,s;k,s) = Mi(k,s;k,-s) = 0. (E2) 

For the out-of-phase part, £ - (q) is linear with q, 

M Q (k,s;k,s) = M (k,s;k,-s) = 0. (E3) 

Therefore, all the zero-order terms vanish. 

First-order: We write the matrix elements of mo- 
mentum scattering using the basis states of the X point. 
Here, the effect of spin-orbit coupling can be safely ne- 
glected. As a result, in the expansion of |k, s) [Eqs. ([26)) - 
(|27p] we omit B(k) and keep only the coefficients of A(k). 
These coefficients and the translational part are then lin- 
earized around K = (ki + k 2 )/2 providing 



A(ki 



cxp(ik' 12 -r) 



A(k ) T q/2- V fc A(k) +.. 

o,i > ±2£l±. p *' 



2E„ 



~2E„ 



g,X jL, - LJ g,X_ 

exp(zk[j • r)(l =p iq ■ r/2 



, (E4) 
, (E5) 



where K was replaced by the wavevector of the valley 
center (ko). The error brought into the matrix element 
with this replacement is of higher order (quadratic) in q. 
We substitute these linearized forms into Eq. (|2"o]) and the 
resulting states are then plugged into the in-phase and 
out-of-phase part of Eq. (JSJ). This procedure identifies 
three linear terms with relatively large coefficients, and 
the matrix element up to leading order reads 

3 

M(kr, s; k 2 , s) = X ( n } + M« = E 7 mJ 

i m ,i = $-(a)-(xl\vv-\xl) 
p 



Im.2 — 



2E„ 



E [^ + (q)-(^l|V^ + l^i 1 >-cc.] 



£=x,y 



Im,3 = -| + (q) • (Xl\(iq-r)VV + \Xl) , 



(E6) 



where c.c. denotes complex conjugate. These conceiv- 
ably dominant terms are further examined to see if they 
are kept by selection rules. 

The integral forms of intravalley momentum scattering 
are further restricted by group theory. Both the basis 
states couplings and interactions are identified as parts 
belonging to IRs of G\ 2 . Using Eqs. (fTg )) - (|20b[) and the 
discussion that follows, we write the following decompo- 
sitions, 



(x\\...\x\) 
{x{\...\x\) 

VV+ 

w_ 



X 1 ®X 1 = Mi © M 3 ©M*.© (E7a) 

X 4 ® Xi = M 5 © M' b , (E7b) 

(M3©M 5 )<8)Mi =M3©M 5 , (E7c) 

(M 3 © M 5 ) $ M 2 = M 4 © M' 5 . (E7d) 



In writing the first line, we have used Eq. (|19a[) to cross- 
out IRs in which only off-diagonal coupling between ba- 
sis states is possible (i.e, between Xf and Xj). The 
first integral, I m j in Eq. (|E6|) . couples X\ states via 
a vector-type operation on the antisymmetrical poten- 
tial part (VV_). Since the respective decompositions in 
Eqs. (|E7a[) and (|E7d[) have no mutual IR, this type of 
coupling vanishes. We get 



In 







(E8) 



The second integral, 7 mj2 in Eq. (|E6I) . couples X\ and X 4 
states via a vector-type operation on the symmetrical po- 
tential part (VV+). M 5 appears in the respective decom- 
positions [Eqs. (|E7b[) and (|E7c[) ] . Following the discus- 
sion that precedes Eq. (|21[) we can find the invariant inte- 
grand form. Then, by noting that V+ can be replaced by 
the crystal potential (set ki i2 = (0, 0, l)27r/a in Ref. [93l|). 
the integral is analytically solved using Eq. (|E1[) . 



Ira. 2 — 



im P 2 



[&£(q) + 9y£y (q)] 



(E9) 



The third integral, I m s in Eq. (IE6[) . couples X\ states 
via a second-rank tensor [r<g>VV+(r)]. To write a de- 
composition expression, it should first be casted into a 
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symmetrized form, 
r©W+ = ^(r - T«)®VKt(r - r a ) + r® W_ . (E10) 

As shown before, VV_(r) cannot couple between Xj 
states. The symmetrized sum term on the right-hand 
side transforms as 

(M 3 © M 5 )® 2 = (M B <g> M 5 ) © (M 3 © M 3 ) © 2(M 3 © M 5 ) 
= (Mx © M 2 © M 3 © M 4 ) © (Mx) © 2(Mg). (Ell) 



Combined with Eq. (|E7ap . we see that the two Mi 
in Eq. (|ETT|) contribute to 7 m , 3 in Eq. (|E6j). There 
should be two independent parameters associated with 
the two Mi's. Out of the nine tensor components of 
S a ( ri — Ta^dV^tij — T a)/drj, one independent param- 
eter originates from the product of longitudinal compo- 
nents (i = j = z) and belongs to M 3 © M 3 = Mi, and 
the other from the sum of transverse component products 
(i = j = x plus i = j = y) and belongs to the Mi out of 
M$ © M5. Putting all of the pieces together, the overall 
matrix element of intravallcy momentum scattering reads 



M(ki,s;k 2 ,s) 
irrioP 2 



(E12) 



h 2 

x,y,z A,B 

; E<^ + E 



(q) + (q)] 



dV at {v- 



<9r. 



Deformation potential theory provides a concise ap- 
pearance for the matrix element of intravalley momen- 
tum scattering , 101 ' 107 ' 129 



E< 



k Vjj k„ < ,,.:q 



(E13) 



where deformation potential operators and strain tensor 
(e) elements at the long-wavelength regime are 



T> jk 
ejfc(q) 



PjPk 



lim 

g->o 



dV e [(l + e)-r] 



mo e-m dejk 
ife^(q) + %e/(q))/2. 



(E14) 
(E15) 



V e is the crystal potential under strain e. Here ejfc(q) 
in Eq. (|E15[) is stripped out of the amplitude factor 
y/H/[2pV(jj(q)]y/n(q) + 1/2 ± 1/2 in order to compare it 
with the expression of Eq. flSJ). To relate with Eq. (|E12[) 
we substitute Eq. (|E14|) and (|E15j) into Eq. (|E13|) and 
use ^k ( r ) — e' lk ° z ipxl ( r )- The second term in leads 
to exactly the last term of Eq. (|E12[) with an opposite 
sign. For the kinetic term in Djj, we write 



(A 1 \p 2 \X\) = {x\ \ Pj \x n ) (X n \ Pj \x±) . 



mlP 2 



, (E16) 



where j = {x, y} considering the dominant coupling of 
Px,y between X\ and valence X4 states*^ p 2 does not 



have this dominant coupling. Therefore, the first term 
in Vjj is exactly the second line of Eq. (|E12[) with a 
minus sign. Physical results do not change upon this 
global minus sign. Of the diagonal deformation potential 
constants Djj = ( ko [Djj | ko), there are two independent 
values D zz and D xx = D yy (in z valley). Dilation and 
uniaxial deformation potentials are related to them via 
= D xx and a u = D zz — D xx . 

All in all, we have shown the equivalency of our pro- 
cedure with known deformation potential theory for in- 
travalley momentum scattering. 



Appendix F: some details in intravalley & g-process 
spin flip matrix elements 

1. intravalley 

We invoke group theory and evaluate the integrals 
in Eqs. (|42b[) - (|42e[) . The out-of-phase matrix element, 
Euv^i,^! includes a single dominant Elliott process. 
Its coefficient product, irj/Ac, comes from A x i and 
dB x2 i/dk x ( y y Its integral reads (Xl\ W- \ X\ ). Using 
appropriate decompositions 

(Xl\...\Xf) ->X 1 (g>X 1 =M^ffiM 3 . / © M 4 © M' 2 (Fla) 
VV- -> (M 3 ffiM s )®M^ = M A ®M' 5 , (Fib) 

we see that these states can be coupled by the longitu- 
dinal component of VF- (transforms as M4). We have 
crossed out IRs in which only diagonal coupling is possi- 
ble (e.g, between X\ and A 1 ). All in all, one combination 
is kept 



dz 



XjVife +<&,)£" (q).(F2) 



The remaining non-vanishing matrix elements relate 
to the in-phase potential. The sum ^ v I^^a in- 
cludes a single dominant Elliott process. Its coefficient 
product, —iPr)/E g X Ac, comes from dB x2 > /dk x y and 

dA x *(v) /dk x ( y y Its integral reads (Xf \W + \X%' V ). Us- 
ing appropriate decompositions 

(X? I ...\X{) -> X 4 © X x = M 5 © M' b , (F3a) 
VV+ -> (M 3 © M B ) © Mi = M 3 © M 5 , (F3b) 

we see that these states can be coupled by the transverse 
components of VF+ (transform as M5). Similar to the 
derivation of the second momentum integral [Eq. (|E9I) ]. 
we reach at an analytical result for the third spin integral, 



irjP 2 rao 
2A c h 2 



(fe-^)[^+(q) + 9^(q)]-(F4) 



The sum v I^^A also includes a single dominant 
Elliott process. Its coefficient product, —irj/Ac, comes 
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from A x i and dB x2 > /dk x ^ y . Its integral reads (Xf |r® 

VV+lXj 1 ). The interaction is first casted into Eq. (|E10j) 
where both parts (tensor and antisymmetric potential) 
can couple between the states. Each of the respective 
decompositions in Eqs. (|E11|) and (|Flbj) share a com- 
mon M4 IR with Eq. (jFlap . Using the transformation 
properties of M4, the resulting integral reads 



iv(Vx - iq y ) 



^c~ V ^ 



2' 



dV- 



dz 



q-rC(q) (F5) 



The sum ^ ^ i^^s includes dominant Elliott pro- 
cesses coupled by VV+. One of the dominant products 
comes from A x i and 2d 2 B x2 > /dk z dk x ( y y The corre- 
sponding integral (X 2 \W+\Xl) vanishes, for that there 
are no common IRs between Eqs. JFTJ) and (|F3bJ) . Other 
dominant products come from Avi and d 2 B VI ( V ) Idk,,,^, 

and from B x<c ( v ) and 2dA x ^jdk x dk y . Repeating the 
analysis that led to Eqs. (|E9|) and (|F4|) . we reach at 



ij.i' 



—iP 2 m 
2A c h 2 



-W-rj)q x q y fe + (q) - *£+(q)] } ■ (F6) 



The only Yafet process is included in I ^^,5 with coef- 
ficient product 4P 2 / (E g x Ec) that comes from A x i and 
2d 2 A x2 , /dk x dk y . Its integral reads (XflVV^jXj 1 ). 
The fc-independent part of V|° k docs not contribute 

to this matrix element: (X 2 |V(VV+ x p)\Xl) van- 
ishes by time reversal symmetry [Eqs. ()22b|) and (|22d|) ]. 
The relatively small interaction V(VU+ x ftko) can cou- 
ple Xl with X\ states. The interaction transforms as 
second-rank tensor and belongs to (M3 M5)® 2 . From 
Eqs. (|E11[) and (jFlap . we see that M4 is the common IR. 
To be specific, it is the component d 2 V + /dxdy that be- 
longs to M 4 . In spite of the relatively small magnitude, 
we still give its expression explicitly, 



c 



AttIqC 2 



Xf 



d 2 V + 



dxdy 



Xl), (F7) 



and shall find it a compensating part in leading to a 
concise result. 

Given these dominant contributions for spin relax- 
ation, we may attempt to relate them to some known 
physical quantities. First of all, we find out that the 
terms in Eq. (|F4|) and Eq. (|F6|) (excluding small if — r) 
part) can be combined and they share a common factor 
iv(Qx — i<ly)(qx£y + Qy£x)/ 'Ac with the second term in 
Eq. (|F5[) . Secondly, the left-out rf — r\ part in Eq. (|F6[) 



is found to be exactly compensated by the Yafet part 
[Eq. (|F7]l]. It is shown with the help ofi2i 



4?7JqC 2 



xl 



d 2 V + 



dxdy 



xl 



2Pm a 
H 2 ' 



(F8) 



At this phase, the total of Eqs. (|F4[) - (|F7[ ) closely resemble 
to a deformation potential constant D' = (X 2 \D xy \Xl) 
defined in Eq. (|E14[) . up to a constant prefactor. How- 
ever, care should be used due to its off-diagonal na- 
ture. Contrary to diagonal deformation potential that 
appeared in the momentum scattering (e.g., dilation and 
uniaxial), 9V^[(1 + e) • r]/de xy contains a part induced 
by internal displacement #i2I Thus, the sum of Eqs. (|F4|) - 
(|F7[) alone is not sufficient to form a complete deforma- 
tion potential constant. Out-of-phase phonon polariza- 
tion vector £ , on the other hand, can be expressed in 
terms of £ + and internal displacement for small q (see 
Table HT1 and its discussion). Replacing £~ in Eq. (|F2|) 
with equivalent £ + terms (from the fourth row in Ta- 
ble HI)) , we obtain the final intravalley spin- flip matrix 
clement expression in Eq. (|45|) with Eq. (|46l) of the main 
text. 

Finally, we mention an alternative choice to derive 
these results with a basis states of the A axis at ko. 
In this case, third-order perturbation theory expanded 
around k = ko gives a similar result (with the help 
of wavevector-order analysis and the appropriate adjust- 
ments to the space group of A-axis). From our results, we 
are able to conclude that the leading contributing term 
with this alternative approach is 



^(A 1 ||q.p|AU)(Am^ |A 2 ,t> 
r q-p|A|t)(A|tll + (q)W 



m Eg^ Ac 
x(A 2 M 



+e xy {q)V xy A^), 



where i,j — {x,y}. p, |A) and E g ^ denote, respectively, 
the momentum operator, the spin-independent state and 
the energy gap between conduction and valence bands 
at ko. £ X y(si)Vxy, with Vjk defined by the last term in 
Eq. (|E14[) . is a partial combination of in-phase and out- 
of-part interactions. 



2. g-process 

We analyze the terms I E and I Y in Eq. and de- 
rive their expressions in detail. Following the reasoning 
that led to observation 2 of the intravalley case, the El- 
liott part (I ) has a dominant coefficient ir\j Ac from 
dB X 2' /dk XiV and A x i. However, the basis states associ- 
ated with this coefficients are (Xj| and \Xl), as can be 
inferred by the index arrangement of // [below Eq. (|53[) ]. 
The Yafet part (I Y ) has a dominant coefficient —P/E giX 
from dA x x(y) /dk x ( y ) and A x i between (X%^\ and \Xl). 



28 



In spite of this coefficient disobeying observation 2, the 
Yafet part is kept and will be shown not too small com- 
pared with the Elliott part. 

Having identified the important coupling between ba- 
sis states, we can determine which of the phonon modes 
dominate the <?-process spin relaxation. It is deter- 
mined by applying selection rules connecting opposite 
points of the A star. The Elliott part is analyzed first. 
From Eq. (|53[) . a first Brillouin zone phonon wavevec- 



tor q = —2kg rj (0, 0, 0.3)27r/a is needed to conserve 



the crystal momentum. With tj)* xl = ip X 2' , the char- 
acter of the state product has the following identity, 

X%X% = xfP = (Xlf°y = (xlT'r- The last equal- 
ity is obvious if one refers to the character table of the 
A group (e.g., Table V of Ref. [H|). Thus, the Elliott 

— 2k' 

part involves a phonon with a character % A °. IR of Ai 
with a first Brillouin zone wavevector is identified with 
the LA mode. Similarly, the coupling of the Yafet part is 

— 2k' 

between V'Ai.fco and ipA 5 ,k , and it leads to Xa ° f° r the 
interaction. The spin-orbit potential transforms as A5, 
and XAi(a^)Xa 5 = XA 5 - Thus, the Yafet part involves a 
phonon mode that transforms as Ai or A' 2 , which corre- 
sponds to LA or LO mode, respectively. The LO phonon 
has much larger energy than that of the LA phonon and 
we drop it from the leading Yafet contribution. In conclu- 
sion, the LA mode dominates both the Elliott and Yafet 
coupling. 

The leading-order matrix clement of the g-process 
[Eq. (|53"j) ] can be related to some form of a deformation 
potential parameter. It can be done since the coupling 
can be expressed between basis states of the same X 
point. The Elliott part is explicitly written as 



-2r? 
A^ 



2ik'z ( e + 



LA,z 



8V+ 

~dz~ 



LA ' z lh 



(F9) 



which belongs to 



M'^ 2 ® M3 © M3 © M3 ® M3 ® M, 



As shown by Eq. (|E7a|) , operators that couple X\ states 
belong cither to Mi or M3. Operators that are all even 
under time reversal [Eq. (|22b[) ] are further restricted to 
Mi for even-parity states [Eqs. (|19aj) and (|22d|) ]. Thus, 
the allowed coupling interaction is J2 a ( z ~ T a,z)dV a t(r — 
T a )/dz. Similarly, we find that there is no comparable 
out-of-phase contribution. The resulting Elliott part of 
the matrix element reads 



D zz {K x -iK y ), 



(F10) 



where we have used £laz( — 2kg) ~ 1 and D zz in 
Eq. dEil) . 

For Yafet part, the unity leading term in the ex- 
pansion of e 2lk o z results in a nonvanishing integral. It 
is the reason that we have kept this part in spite of 
disobeying observation 2. The symmetry properties 
of d{VV+ x p) x ( y )/dz follow M^® M' b = M 5 which 
also appears in the decomposition of the states product 



[Xf x) ®X{- see Eq. 
coefficient of the state expansion, 
element of the Yafet part reads 



c])]. Together with the leading 



-P/E g 



I Y = 



2P 



D so (K x - iK y ), 



the matrix 



(Fll) 



where the scattering constant integral is denoted by 



where we have used the approximation that nas °nly 
a nonzero longitudinal component (it is exact when the 
wavevector is on the A axis). To find the nonvanishing 
part of this integral, the integrand is converted into parts 
belonging to IRs of the A-point space group. The new 
feature that emerges in a g-process is the phase factor 
e 2zk' z = j + 2ik' oZ - 2k' 2 z 2 + O((k' z) 3 ). For the in-phase 
part, dV+/dz and zdV+/dz belong to M' z and Mi © M 4 . 
The higher-order term is 



z2 dv+ 

dz 
+2t z 



2 dV at (r-T a ) 2 dV+ 

z - T a,z) ^ h T. 



,dV a t{T-T A ) , , 
\Z-T A ,z) 5 (z-T B ,z) 

oz 



dz 

dV at (r~T B ) 

dz 



4c 2 m,Q 



XI 



d(vv+ x P )„ 



dz 



6.7 meV • 2?r/a 



(F12) 



with its value calculated from EPM. The out-of- 
phase part leads to a slightly smaller coupling inte- 
gral («4 meV • 27r/a), but with a small polarization 
^£" Az (— 2fc ) m 0.2. This property renders its contribu- 
tion too small compared with the leading Elliott part. It 
is therefore neglected. 

All together, the leading g-process matrix element is 
Eq. ([54]) with Eq. flS5j in the main text. 
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do not depend on the specific choice of orthogonal polar- 
izations. 
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ther bands explicitly into consideration, the kinetic term 
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if our state basis includes these further bands. 
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